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ABSTRACT 

The primary objectives of this paper are 
pedagogical: to provide a reliable semi-subjective transformation 
procedure that might be used without, difficulty by beginning students 
in factor analysis; to clarify and extend the existing knowledge of 
oblique transformations in general; and to provide a brief but 
meaningful explication of the general obliguimax. Implicit in these 
three objectives is the fourth objective of presenting a paper that 
might be profitable for both the beginning student and the factor 
analytic theoretician- The first section, primarily background, 
discusses one of Thurstone's methods of determining oblique 
transformations. The second section discusses certain theoretical 
aspects of the general obliguimax to provide a basis for the 
development and understanding of the simplified obliguimax. Finally, 
the semi-subjective simplified obliguimax transformation is developed 
and discussed within the context of Thurstone's box problem. 
(Author/AE) 
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Once a factor analyst has determined a factor space for a complex of 
variables he usually desires a basis for interpreting the factor space. 
Traditionally the basis for interpreting a factor space has been made 
relative to either orthogonal (uneorrelnted) factor or oblique 

(correlated) factor axes. That is to say, the "loadings’ of the initial 
factor loading matrix are transformed through the use of either an 
orthogonal or oblique transformation procedure. 

Any oblique transformation solution encompasses three matrices. The 
general names given to these three matrices are the pattern matrix, the 
structure matrix and the factor intercorrelation matrix. The entries of 
a structure matrix represent the perpendicular projections of the variable 
vectors onto the oblique factor axes and with an appropriate scaling are by 
rows the correlations of' the variables with the factors. The entries of a 
pattern matrix represent the parallel projections o£ the variable veelqis 
onto the oblique factor axes and are by rows the stand irdix id re wv ns ion 
weights of a regression equation describing each of tno observed variables 
in terms of the correlated factors. The entries of the factor 
intercorrelation matrix are just the cor re let ions between the factors. 
Although most papers on oblique transformations do not deal specifically 
with orthogonal factor axes, one may regard an orthogonal transformation 
solution within the more general framework of oblique solutions. In the 
general oblique framework the orthogonal tr uisformati'M solution mi.:ht bo 
thought of as a special solution in which- the factor irta re rr . lotion 
matrix is an identity matrix. Wien the factors are un cor related the 
parallel projections and perpendicular projections of the variable vectors 
are identical, thereby resulting in a structure matrix that is identical 
to a pattern matrix. 
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There have been in the past two schools of thought pertaining to 
the type of interpretation that should be made for an oblique solution. 
One mode of thought was based upon the work of Louis Thurs tone ( 1.947) 
while the other was based upon the work of Karl Helzinger (ho 1 zii vr 
and Harman, 1941). The Thurs tone-Ilolzinger differences stem, 
philosophically from the idea of invariance and geometrically from the 
definitions of factor axes used in obtaining the final solution 
matrices, which implicitly determine whether it in the pattern or 
structure matrix that should be used for the final interpretation. The 
Thurstone (1947) school of thought bases the final interpretation of 
the oblique transformation on a structure matrix while the Ilolzinger 
school (Holzinger and Harman, 1941) uses a pattern matrix for the final 
interpretation. 

Ilolzinger defined his solutions using primary vectors, those 
vectors formed by the intersection of hyperplanes. Thus, Ilolzinger* s 
solution matrices for interpretation were the primary pattern and the 
primary intercorrelation matrices. The loadings of the primary pattern 
matrix are defined geometrically as the parallel projections of the 
variable vectors onto the unit length primary vectors. 

Thurstone (1947) defined his solutions using reference vectors, 
those vectors defined as normals to hyperplanes. Thurstone was 
concerned with the perpendicular projections of the \virlable vectors 
onto the unit length reference vectors. Although Thurstone was 
interested in the reference structure matrix, it is interesting to note 
that he usually reported the primary intercorrelation matrix along with 
the reference intercorrelation matrix. 
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Typically a solution matrix is desired which will have scientific 
meaning and interpretability. Scientific meaning and interpretability 
are facilitated when some of the entries of the solution matrix are very 
high and the remaining entries are zero or near zero, A zero entry in a 
solution matrix may be thought of as a vanishing projection, thus the 
objective in obtaining a solution facilitating scientific intern rotation 
is one of maximizing the number of vanishing projections. (This concept 
is frequently referred to as "simple structure", however this term is 
somewhat misleading and will not be used in this paper.) Either the 
number of vanishing perpendicular projections or the number of vanishing 
parallel projections must be maximized, inasmuch as both types of 
projections cannot generally be maximized within the context of a single 
(either primary or reference) system. That is to say, the zero 
"loadings" in the pattern and the structure matrix cannot both be 
maximized within a single system. The matrix in which the vanishing 
projections are to be maximized is dependent upon the interpretation 
that one wishes to make from the final solution. If the interpretation 
is to be made in terms of the correlations between the variables and 
the factors then the vanishing perpendicular projections of the 
variable vectors onto the unit reference vectors, the near-zero 
entries of the reference structure matrix, should be maximized. If one 
wishes to treat the observed variables as dependent variables and the 

A 

factors as independent variables, then the vanishim r . parallel 
projections of the variable vectors onto the unit primary vectors, the 
near-zero entries of the primary pattern matrix, should be maximized. 

In the past, with several exceptions, most attempts to develop 
analytic oblique procedures have followed the Thurstonian mode of 
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thought, maximizing the vanishing perpendicular projections of the 
variable vectors onto the unit reference vectors. The reason for this 
is not particularly clear, but either the Thurstonian anpronch is l-'ss 
complex than the HoLzinger approach or factor analysts find the 
reference structure matrix an easier matrix to interpret, in keeping 
with tradition this paper will follow the Thurstonian model, howew r it 
should beco. : . \. .rent to the reader, as a result of reading this pa;v.r, 

that the c\s v .* 1 ' model in this paper was somewhat arbitrary as both 
types of solur* . , may be computed with ease using the now transformation 
procedure developed and presented herein. 

The object! .e of this paper is to acquaint the reader with certain 
aspects of the methodology, properties and nature of the general 
obliquinax transformation (Hofmann, 1971). ”or pedagogical and illustrative 
purposes one of Thurstone’s methods of determining oblique transformations 
is presented and then modified to produce a simplified ver/i/m of the 



obliquimax which is referred to as the simplified obi iquimax. 

The simplified obliquimax is unique in the sense that it is a 
semi-subjective transformation procedure that depends neither on an oblique 
analytic simple' structure criterion nor graphical techniques to 
determine the oblique transformation solution. It provides a 
conceptually simple yet reliable oblique transformation proc-nlure for 
most sets of data. However, it is not in ; : ; to be pr.ict !«mJ. 
working model. It is the general obliquimax that is the p>* 'tical 
model. 

This paper is composed of three sections. In the first section, 
(Section I), Thurstone's (1947) method of determining subjective oblique 
transformations is discussed within the context of two-dimensional 
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sections. Although nothing in the first section is new, it defines and 
describes geometrically the matrices and terminology traditionally used 
in the Thurstonian type oblique transformations. The first section also 
establishes an algebraic model for determining subjective oblique 
transformation solutions through the use of an iterative method. 

In the second section, (Section II), the general obliquimax is 
briefly discussed with respect to the matrix equations defining an 
oblique solution. Special emphasis is placed on the matrices of 
direction numbers and solution matrices expressed within the metric of 
the original factor solution. Several important similarities between 
the Thurstone and Holzinger solutions are noted. This section does not 
provide a detailed discussion of the general obliquimax inasmuch as it 
is included only to provide a basic theoretical rationale for the 
development of the simplified obliquimax in the third section of the 
paper. 

In the third and final section, (Section III), of this paper the 
simplified obliquimax is presented. Thurstone' s initial matrix of 
direction numbers is defined as a symmetric matrix a priori without the 
use of planar plots. All subsequent iterative stage solution matrices 
are expressed within the metric of the original factor solution and 
defined in terms of an orthogonal transformation of the original 
factor solution and exponential powers of the initial symmetric matrix 
of direction numbers. Conjectures are made about certain now aspects 
of the geometry of oblique solutions within the framework of the 
direction numbers of the simplified obliquimax. 

In Sections I and III a set of illustrative data is used to 
'clarify the discussion. Iterative solutions for this data set are 
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determined by Thurs tone’s method in Section I and by the simplified 
obliquimax in Section III. 

Aside from acquainting the reader with the general obliquimax, this 
paper should clarify and extend certain theoretical aspects of oblique 
solutions in general. 



Section I 

Thurs tone's Method Of Determining Oblique Transformation 
Solutions By Two Dimensional Sections : 

As previously mentioned Thurstone sought to define an oblique 
solution with respect to perpendicular projections onto the reference 
vectors, thus his solution matrix of interest was the reference 
structure matrix. Thurstone (1947) presented several methods of 
oblique transformation: plotting the normalized variable vectors onto a 

hyper-sphere, two dimensional sections and by three dimensional sections. 
His first method was quite subjective while his other two methods were 
primarily analytic and algebraically the principle involved in both 
methods is the same. In this section his algebraic principles will be 
used and discussed within the context of two-dimensional sections. 
(Although there are numerous modifications and rewordings this section 
is taken directly from Thurstone (1947, p. 194-224). Reference to 
Thurstone (1947) is implicit throughout this section.) 

Assume some factor loading matrix, F, defining the perpendicular 
projections of n variable vectors onto v mutually orthogonal factor 
axes. The r axes are arbitrarily orthogonal axes as determined by the 
initial factoring method. For illustrative purposes Thurstone’s 
technique will be discussed within the framework of his classic box 
problem (Thurstone, 1947, pp. 140-144). The centroid solution for the 
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box problem is reported in Table 1. For this particular set of data 
(n = 20) and (r = 3). 



Table 1 



Centroid Solution* of Box Problem 
Matrix F 







c o 






h 


B i 


c i 


1 


.659 


-.736 


.138 


2 


.725 


.180 


-. 656 


3 


. 665 


.537 


.500 


4 


.869 


-.209 


-.443 


5 


.834 


.182 


.508 


6 


.836 


.519 


.152 


7 


.856 


-.452 


-.269 


8 


.848 


-.426 


.320 


9 


.861 


.416 


-.299 


10 


. 880 


-.341 


-.354 


11 


.889 


-.147 


.436 


12 


' .875 


.485 


-.093 


13 


.667 


-.725 


.109 


14 


.717 


.246 


-.619 


15 


.634 


.501 


.522 


16 


.936 


.257 


.165 


17 


.96b 


-.239 


-.083 


18 


.625 


-.720 


. 106 


19 


.702 


.112 


r r* r*. 

-. oju 


20 


.664 


.536 


. 483 



o 

ERIC 



*Thurstone, 1947, p. 194. 

The r axes may be denoted as A , £ , and C . These arbitrary axes 
are regarded as fixed in position. The problem is to select by u 
successive approximations the unit reference vectors. /•'. , and 6\ . 
such that the number of variable vectors with zero perpendicular 
projections onto these unit reference vectors is a maximum. 

In starting the transformation procedure the r unit reference 
vectors are assumed to be mutually orthogonal and collinear with the r 
arbitrary orthogonal axes. That is to say, Aj is orthogonal to Bj and 
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Cj and it is collinear with A Q . This may be observed in Figures 1, 2 and 
3 which represent the planar plottings of the variable points with 
respect to. the initial factor axes as reported in Table 1 (Figure 1 being 
the points as defined by the first two columns of Table 1). Note that 
A 0 and A ^ are the same axis. 

Figure 1 

Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
A 2 and 2^, Projected Onto Plane A.Pj 







The locations of the reference vectors may be defined with respect 
to the fixed orthogonal frame through the use of the matrix of 
direction cosines V . The subscript of V • refers to the th given 
positions of the reference vectors (the iterative stage). When v is 
unity the reference vectors are collinear with the fixed orthogonal 
frame and Vj is an identity matrix. The columns of Vj give the 
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direction cosines of the initial locations (1) of the unit reference 
vectors with respect to the fixed orthogonal frame.* For any matrix of 
direction cosines, V u , the entry • refers to the cosine of the angle 
of inclination between the unit reference vector j and the original 
fixed axis i. 

The points in Figure 1 show the configuration of variable vector 
termini as they would appear when projected orthogonally onto the plane 
of A If vector A ^ is transformed in the plane AjBj to the position 
of A fa, it will determine a plane (of hyperplane if r > 3) which will 
intersect the plane AJB^ in the line which is marked B^-primary. The 
vectors associated with 1, 13 and 18 will have near-zero projections 
(vanishing projections) on the vector A!>. It is important to note that 
the B^-primary passes through the group of points 1, 13 and 13. 

Similarly the given position of B^ can be transformed in the A^B^ plane 
to B| and its associated plane will intersect the plane A 7 B^ in the 
line marked A^-primary. The ^-primary passes through the group of 
points 6, 9 and 12 and their variable vectors have vanishing projections 
on Bo, 

In transforming A ^ and Bj to the positions A L and Bo respectively 
new positions have been estimated graphically for the reference vectors 
A and B such that the number of variable vectors having vanishing 
projections has increased. • It is important to note here that /l- and B^ 
are in part bases and altitudes of right triangles whose hypotenuses 
are A^, B|, A^-primary and B^-primary. The vectors A£, Bn, AL- primary 

*The initial factor loading matrix F is assumed to represent the 
first iterative stage of Thurstone’s solution. Technically F should be 
subscripted as Bp however for convenience the subscript 1 is omitted. 
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and ^-primary are not of unit length. The prime is used to signify that 
AL and Bl are '''long reference vectors" and that Al- -primary and ZJ£-primarv 

Cl Ct tmj 

are "long primary vectors". 

The geometric discussion in this section is not quite the same as 
that presented by Thurstone. It is hoped that by including the long 
primary vectors as well as the long reference vectors some of the 
geometric similarities between the Thurstone and Folzinger solutions 
will become evident. Technically a reference vector is orthogonal to a 
hyperplane of (r - 1) dimensions. In any r-dimensional space there are 
r hyperplanes of ( r ~ 1) dimensions, and therefore r reference vectors. 
The intersection of (r - 1) hyperplanes defines a primary vector, 
therefore in any r-dimensionai space there are v primary vectors. The 
vectors drawn orthogonal to a hyperplane will necessarily be orthogonal 
to any vectors contained within the hyperplane., thereby implying that 
each primary vector must be orthogonal to (r - 1) reference vectors. 
Therefore each primary vector is correlated with only one reference 
vector. Orthogonal to the one hyperplane not containing the primary 
vector is that reference vector. Within the context of this paper each 
reference vector is referred to by a subscripted Roman letter. The 
Roman letter may be thought of as representing the hyperplane to which 
the reference vector is orthogonal. Each long primary vector is 
referred to by a subscripted Roman letter. The Roman letter associated 
with the long primary vector may be thought of as representing the 
hyperplane which does not contain the primary vector. Thus for the 
illustrative example long primary is orthogonal to all long reference 
vectors with the exception of long reference vector A!). This discussion 
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may be generalized to any number of dimensions* and will be presented 
algebraically in Sections II and III. At this point we only wish to . 
call the reader's attention to the long primaries in the plotted figures 
and to note that they are the Ilolzinger (Ilolzinger and Harman, 1941) 
long primaries. 

The coordinates of the termini of An, S', A 4 -primary and # '-primary 

^ Ci 

can be defined with respect to the fixed orthogonal axes A 0 , B Q and C 
or with respect to A^, Bj and Cj. Only the coordinates of the long 
reference vectors will be discussed in this section. The termini of the 
long reference vectors AJ and B ^ are linear combinations of Aj and B 
Specifically:- 
A£ = l.OOA^ + .90^; 

Bij » . 50A^ - 1.00 B r 

The coordinates of the terminus of AJ with respect to AL and B ^ are 
(1.00, .90) and the coordinates of B ^ are (.50, -1.00). The use of A Q 
and A^ may be somewhat perplexing to the factor analyst unfamiliar with 
Thurstone's methodology. For subsequent iterations the role of A 0 as 
opposed to the role of the previous position of the reference vector, 

A u _2 which is A'j for the first iterative stage, will become much clearer. 

In Figure 2 the first and third columns of F have been plotted. 

The vector C n has been transformed to CL such that A '-primary passes 
through the group of points 8, 11 and 18. The variable vectors 
associated with these points will have vanishing projection.; on the long 
reference vector C '. The coordinates of the terminus of CL with respect 

a ^ 

to A 2 and are (.40, -1.00). 

C L = . 40A - l.OOtf 

u 1 1 

*When r < 4 the hyperplanes become planes. 




12 



12 



* 



i 



i 



» 



Figure 2 

Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
A j and Cj> Projected Onto Plane A jC 1 




Note that A is not transformed in the AjC ^ plane. The axis A j in the 
A j plane represents the projection of long reference vector onto 
the plane. 

In Figure 3 the second and third columns of F have been plotted. 

In the B^Cj P lane the axes ^ 2 an( * ^1 are n0t trans f° rine< *> thus, in this 
plane B^ becomes long. reference vector B^ and becomes long reference 
vector The axes B^ and represent the projections of long 

reference vectors B^ and respectively, onto the B-C? plane. 

Although Thurstone has described the transformation process in 
two-dimensional sections the long reference vectors are theoretically 
linear combinations of all three of the initial reference vectors. They 
are also a linear combination of all three of the original axes. 
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Figure 3 



Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
B 1 and C 1t Projected Onto Plane b ,C « 
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The equations for A^ t B 2 and are correctly expressed as: 

* 1 . 00 ^ + . 90 ^ + 0.0 C 1 ; 

B' 2 = . 50,4 ^ - l.OOB^ + 0.0 C 2 ; 

C f 2 = .404^ + 0.0 B 2 - 1.00C- . 

(Thurstone, 1947, pp. 197-198) 

The coordinates of the termini of the new long reference vectors. 
(2) with respect to the previous unit length reference vectors (.7) are 
referred to as .the direction numbers of the reference vectors with 
respect to the previous reference vectors. The matrix of such direction 
numbers will be referred to as S , where the subscript u refers to the 
reference vectors of the present iterative stage and the subscript m 
refers to the previous iterative stage (m - u - 1) . The matrix of such 



direction numbers for the second iterative stage of the illustrative 
example is reported in Table 2. 



Table 2 



Direction Numbers of Second Iterative Stage 
Reference Vectors With Respect To Previous 
Iterative Reference Vectors* - Matrix 



(? 





A 2 


B 2 






1.000 


0.500 


0.400 


B i 


0.900 


- 1.000 


0.000 


c i 


0.000 


0.000 


- 1.000 



*It is assumed that F represents the 
first iterative stage. Thurstone, 
1947, p. 198 



The entries of SL by column are the coordinates of the reference 
vectors with respect to the reference vectors of the previous stage. 

A second matrix of direction numbers is used by Thurstone. The 
second matrix, represents by columns the coordinates of the termini 
of the reference vectors of the u - th iterative stage with respect to the 
axes of the fixed orthogonal frame. The matrix L, is the product of the 

iri 

previous matrix of direction cosines, V , post-multiplied by the present 
matrix of direction numbers, S' , of the reference vectors with respect 
to the m - th set of reference vectors. 



L = V S 

u m mu • 

For the second iterative stage V is an identity matrix therefore 

Lg is identical to » however this identity will not hold for 

subsequent iterative stages. 



[ 1 ] 
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As previously noted the primes on the new reference vectors 
distinguish the long reference vectors from unit length reference 
vectors. The length of the long reference vectors can be determined 
through an r-dimensional extension of the Pythagorean theorem: the sum 

of the squared lengths of the legs of a right triangle is equal to the 
squared length of the hypotenuse. Within the present framework the 
length of a long reference vector is determined with respect to the 
fixed orthogonal frame. Each coordinate of the terminus of a reference 
vector is analogous to the length of one of the legs of a right triangle 
whose hypotenuse is the long reference vector. The direction numbers 
of interest for determining the lengths of the long reference vectors • 
are the column entries of L . The squared length of long reference 
vector Ag is just the sum of the squared entries of the first column of 
Lg (For this particular iterative stage Lg = S ^ and the column sums of 
squares may be computed directly from S^* Table 2.). 

Let the non-zero entries of the positive definite diagonal matrix 

2 

D u represent the squared lengths of the long reference vectors determined 

o 

by the w-th iteration. The equation for computing D u is: 

B 2 u - diagonal (L^L U ). [ 2 ] 

For the w-th iterative stage the value <d»* represents the length of 

the i-th long reference vector. For the illustrative example the 

diagonal entries of D 0 are reported in Table 3. Post-multiplying L by 

z it 

-1 

£> will rescale the metric of the direction numbers such that the 

u 

column sums of squares will be unity for the matrix product /.. . 

Within a trigonometric framework the rescaling of the columns of £ is 
tantamount to dividing the length of each leg of a right triangle by its 
hypotenuse, thereby converting the direction numbers to direction cosines. 
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Table 3 



Lengths of the Long Reference Vectors Determined 
in the Second Iterative Stage* - Matrix 



A 

2 


K 

o 




1.345 


0.000 


0.000 


0.000 


1.118 


0.000 


0.000 


0.000 


0.929 



*Thurstone, 1947, p. 198 

The lengths of the long reference vectors have not been changed. The 
metric of the long reference vectors has been changed. In changing the 
metric of the coordinates of the termini of the reference vectors each 



coordinate becomes the cosine of the angle of inclination between a 
particular reference vector and a fixed axis. The rescaling of L by 
D u normalizes the columns of L u to form the matrix of direction cosines 
I 7 . The equation for computing the direction cosines associated with 



the unit reference vectors of the u- th iteration is: 

F = L D ~ 2 . 
u ■ u u 



C3] 



Table 4 

Direction Cosines of Second Iterative 
Stage* - Matrix V % 





A 2 


h 


C 2 


' A o 


.743 


.447 


.371 


B o 


.669 


-.894 


.000 


C o 


.000 


.000 


-.928 



*Thurstone, 1947, p. 198 
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To determine the perpendicular projections of the n variable vectors 
onto the r unit length reference vectors as estimated by the u- th 
iterative 'stage it is necessary to post-multiply F by the matrix V . 

* vi 

Postmultiplying F by V (Equation 4) will transform the initial reference 
vectors into the positions of the u- th estimate of the reference vectors 
and the resulting matrix, F , will be the reference structure matrix for 
the u - th iterative stage. For the illustrative example V is reported 

u 

in Table 4 and F ' is reported in Table 5. 

F = FV S D~ 2 



u 



u 



m mu u 

- 1 



FL D 
u u 



F„, .= FV. 



u ' 



u 



[4] 



Table 5 

Reference Structure Matrix as Determined By 
the Second Iterative Stage* - Matrix 







B 2 


C 2 


1 


-.003 


.953 


.116 


2 


.659 


.163 


.878 


3 


.853 


-.183 


-.217 


4 


.506 


.575 


.734 


5 


.741 


.210 


-.162 


6 


.968 


-.090 


.169 


7 


.334 


.787 


.567 


8 


.345 


.760 


.018 


9 


.918 


.013 


.597 


10 j 


j » 426 


.698 


.655 


11 i 


1 .562 


.529 


-.075 


12 ; 


i .975 


-.042 


.411 


13 


.011 


.946 


.3 46 


14 


.697 


.101 


.840 


15 


;806 


-.164 


-.249 


16 


.867 


.189 


.194 


17 


.558 


.645 


.435 


18 


-.017 


.923 


.078 


19 


.597 


.214 


.864 


20 


.852 


-.182 


-.207 



*Thurstone, 1947, p. 198 
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Once the reference structure matrix has been determined for the u-th 
iterative stage Thurstone (1947, p. 205) suggested that the cosines of 
the angular separations of the new reference vectors be assessed. His 
objective in assessing these cosines was one of being sure that none of 
the reference vectors had coalesced. If two reference vectors have 
coalesced the cosine of their angular separation, their correlation, will 
be unity and it may be assumed that the problem of transforming to 
singularity has occurred. 

The cosine of the angle of inclination between any two reference 
vectors, their correlation, is the scalar product of their paired 
direction cosines with respect to the fixed orthogonal frame. Equation 
5 may be used to determine the matrix of intercorrelations, Y , of the 
reference vectors of the u- th iterative stage. 



J u = 
*u- 



V'V 

u u 

-7 -7 

D X L , L D M 

u u U u 

D‘ 1 SJV , V S D~ 2 

u mu m m mu u 



[5] 



The intercorrelations of the reference vectors as determined in the 



second iterative stage for the illustrative example are reported in 



Table 6. 



Table 6 



Intercorrelations of Unit Length Reference 
Vectors as Determined in the Second 
Iterative Stage* - Matrix Y 0 





A 2 


B 2 


C 2 




1.000 


-0.266 


0.276 




-0.266 


1.000 


0.166 


C 2 


0.276 


0.166 


1.000 



*Thurstone, 1947, p. 198 
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If the magnitudes of the small reference structure values in the 

matrix F u have diminished with respect to F ffJ and if the magnitudes of 

the large reference structure values in I' u have increased with respect 

to F v „ it is reasonable to proceed to the next iterative stage. However, 
rn ..7 

if two unit length reference vectors are highly correlated it is most 
* .prudent to stop at this particular stage. 

Assessing F with respect to F (Tables 5 and 1) and seeing that the 

t u 

magnitudes of the entries in Table 6 are small we may progress to the 
next iterative step for the illustrative example. 

Subsequent Iterative Stages 

Plot a new set of diagrams for all pairs of columns of F u and 

examine them to determine the adjustments which will be made in this 

particular iterative stage. It is important to note here that Thurstone 

plotted F on orthogonal coordinate cross-section paper, even though the 

reference vectors were oblique. Thus, the reference vectors of F u were 

plotted as being orthogonal. The logic behind this procedure although 

basically simple is frequently quite confusing. For conceptual purposes 

it may be assumed that the axes remain invariant and that it is the 

configuration of variable vectors that is being transformed. The 

apparent paradox is associated only with the plotting of the configuration, 

not with the algebra or interpretation of the reference structure loadings. 

(The initial papers on the obliquimax were all written within this 

conceptual framework and a majority of the algebra was also interpreted 

within this framework.) The new sets of diagrams for all pairs of 

columns of F are presented in Figures 4, 5 and 6. 

2 

In Figure 4 the termini of the variable vectors are plotted with 
respect to unit reference vectors A and B . In the third iterative 

a u 

o 
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Figure 4 

Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
A and 3 0 , Projected Onto Plane A r B c 




stage unit reference vector B is transformed to position B' and 

2 3 

A '-primary will pass through the group of variable points 3, 15, 20, 6, 

O 

9 and 12. The variable vectors associated with these points will have 



vanishing projections on long reference vector B'. In the A e 5 0 plane 

v CJ 2 

the new long reference vector B' may be described in terms of A n and B 

3 2 2 

B' 3 = 1.00 B 2 + .10 A 2 



In Figure 5 the termini of the variable vectors have been plotted with 

respect to A ^ and C Unit length reference vector /■ 0 has been transformed 

to A* , Notice that C'-primary still passes through the variable points 13, 
6 3 

1 and 18 and additionally through the points 7, 10, 4, 9, 2 and 14, thus 
increasing from three to nine the number of variables with vanishing 
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Figure 5 

Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
A and C , Projected Onto Plane A 0 C o 




projections on the long reference vector A£. The new long reference 

vector A' in the A C plane may be described in terms of the unit 
3 2 2 

reference vectors A% and 
A' = 1.004 0 - .72 C. 

v 6 6 

In the same figure C has been transformed to C ^ such that the 
A '-primary passes through the group of variable points 3, 5, 15 and 20. 

O 

The variables 13, 1 and 18 retain their vanishing projections on C ’ and 
additionally variables 3, 5, 15 and 20 also have vanishing projections 

on the long reference vector C f . The number of vanishing projections on 

3 

the C reference vector has increased from three to seven in the A^C % 

plane. Long reference vector C’ may be described in terms of the unit 

6 

•reference vectors A n and C . 

2 2 

C' 3 = 1.00 C 2 + .25 A 2 

o 
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Figure 6 

Planar Plots of Variable Vector Termini 
With Respect to Unit Reference Vectors 
B* and C , Projected Onto Plane l\,C 0 

^ u u u 




In Figure 6 the termini of the variable vectors have been plotted 



with respect to unit reference vectors B and C . Unit reference vector 

2 2 

C has been transformed to Cl so that the B '-primary passes through the 
2 6 S 

variable points 1, 13 and 18. The new long reference vector C 1 in the 

3 

B C plane may be described with respect to unit reference vectors B 

a u a 

and C , 

u 

C' = 1.00(7 0 - ..12B„ 

6 a 2 



In the same figure unit reference vector B 0 is transformed to Bl 

a 6 

such that the C'-primary passes through the group of variable points 2, 

6 

14 and 19. The variable vectors 2, 14 and 19 will have vanishing 
projections on the long reference vector B’. The new long reference 



o 
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vector B f in the C B plane may be described with respect to unit reference 

6 a a 

vectors B ^ and C ^ 

B'g = 1 . 00 B 2 - . 20 ^ 

It is not at all uncommon for a reference vector' to be transformed in 
numerous planes during a single iterative stage if (/’ ;• 3), inasmuch as 
there are r(r - l)/2 possible planar transformations. In this iterative 
stage two of the reference vectors were transformed in two pianos, thus 
they must be described with respect to all three unit reference vectors. 

The final descriptions of the new long reference vectors with respect to 
the previous unit length reference vectors A n , B 0 and C are: 

Ci u Ct 

A ' = 1.004. + O.OOfl - .72 C ; 

3 2 ■ 2 2 

B ’ = .10 A + 1.00B - .20 C i 

3 2 2 2 

C' = .25^4. - . 12S +1.00C . 

3 2 2 2 

(Thurstone, 1947, p. 209) 

Thurstone (1947, p. 208) referred to the coefficients of the above 
equations as "corrections'' . The implication was that the new 
transformations of the reference vectors are simply corrections of the 
previous transformation. It is prudent to realize that these corrections 
as reported in -this paper were determined subjectively by Thurstone. He 
has numerous subjective "rules of thumb" that he utilized when determining 
these corrections (See Thurstone, 1947, pp. 207-210; 212-216). 

The coefficients of the three linear equations arc the direction 
numbers for the three new long reference vectors with respect to the 
previous unit length reference vectors. That is to say, the coefficients 
are the entries of the matrix of direction numbers S „, n - 2 and u = 3. 

O 

Cj t/ 

The entries of Table 7 by column represent the coordinates of the 
termini of the new long reference vectors with respect to the previous 

o 
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Table 7 



Direction Numbers of Third Iterative Stage 
Reference Vectors with Respect to the Unit Reference 
Vectors of the Second Iterative Stage* - Matrix 

t4 O 



i 

j 


3 


O 


C' 

3 




t 

1.000 


0.100 


0.250 


*2 


0.000 


1.000 


-0.120 


C 2 


-0.720 


-0.200 


1.000 



*Thurstone, 1947, p. 209 



unit reference vectors. The coordinates of the terminus of B' with 

6 

respect to A B g and C % would be (.10, 1.00, -.20). 

To convert the direction numbers S ^ to the matrix of direction 
numbers with respect to the fixed orthogonal frame , £ Equation 1 is 



used. 



£ 3 = V 2 S 23 



In Table 8 the direction numbers of the new long reference vectors 

A' y B f and O ' , are reported with respect to the axes of the fixed 
3 3 3 

orthogonal framework, A q> B q and C . 

Table 8 



Direction Numbers Of Third Iterative Stage 
Reference Vectors With Respect To The Original 
Fixed Orthogonal Framework* - Matrix L ? 





A 1 

A 3 


B S 


C 3 




.476 


.447 


.503 


B o 


.669 


-.829 


.275 


C 


.668 


.186 


-.928 


o 









*Thurstone, 1947, p. 209 
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The 



entries of Table 8, by column represent the coordinates of 



the termini of the long reference vectors with respect to the axes i , 

B and C. With respect to the original axes /l , £ and C the coordinates 
OO r O’ O 0 

for the terminus of C* are (.503, .275, -.928). 

6 

Table 9 



Lengths Of The Long Reference Vectors Determined 
In The Third Iterative Stage* - Matrix D„ 



A' 

3 


B 1 

3 




1.06 


.00 


.00 


.00 


.96 


.00 


,00. 


.00 


1.09 



As in the previous iterative stage it is necessary to rescale the 

metric of the long reference vectors to that of unit length reference 

2 

vectors. The squared lengths of the new long reference vectors, D , are 

6 

computed through the use of Equation 2, Table 9. Matrix L is column 

U 

normalized to form the matrix of direction cosines, V^ f of the new unit 

length reference vectors with respect to the fixed orthogonal frame, 

0 

Equation 3. The matrix of direction cosines for the third iterative 
stage is reported in Table 10. 

Table 10 




Direction Cosines of Third Iterative 
Stage - Matrix V 

6 





A 3 


S 3 


C 3 




0.450 


0.466 


0.461 




0.632 


-0.864 


0.252 


C o 


0.631 


0.194 


-0.851 



I 



> 



*Thurstone, 1947, p. 209 
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The reference structure matrix associated with the third iterative 
stage, F„ f is computed using Equation 4. The reference structure ma:rix 
F„ is reported in Table 11. 

O 

Table 11 

Reference Structure Matrix As Determined 
By The Third Iterative Stage* - Matrix F ^ 





A 3 


B 3 




1 


-.082 


.970 


.001 


2 


.026 


.055 


.938 


3 


.954 


-.058 


.016 


4 


-.021 


.500 


.723 


5 


.811 


.330 


-.002 


6 


.800 


-.031 


.387 


7 


-.071 


.737 


.510 


8 


.314 


.825 


.011 


9 


.462 


-.017 


.756 


10 


-.043 


.636 


.621 


11 


.582 


.626 


.002 


12 


.642 


-.030 


.605 


13 


-.090 


.958 


.032 


14 


.088 


.001 


.919 


15 


.931 


-.037 


— . 026 


16 


.688 


.246 


.356 


17 


.231 


. 641 


.456 


18 


-.070 


.946 


-.035 


19 


-.024 


.104 


.905 


20 


.945 


-.060 


.026 



*Thurstone, 1947, p. 209 

The intercorrelations of the new unit length reference vectors are 
computed using Equation 5. The matrix of intercorrelations of the unit 
length reference vectors, Y , as determined in the third iterative stage 

O 

is reported in Table 12. 

The magnitudes of the loadings in Table 11 are assessed with respect 
to the loadings in Table 5. Clearly the small loadings are diminishing 
and the large loadings are increasing. None of the off-diagonal entries 
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Table 12 



Intercorrelations of Unit Length Reference Vectors 

As Determined By The Third Iterative Stage* - Matrix Y , 

o 





A 3 


B s 


C 3 


A 3 


1.000 


-0.214 


-0.171 


B 3 


-0.214 


1.000 


-0.169 


C 3 


-0.171 


-0.169 


1.000 



*Thurstone, 1947, p. 209 



of Y are approaching unity. A fourth iterative stage is in order. 

O 

A detailed discussion of the fourth iterative stage will not be 
presented in this paper. As with the previous iterative stages the 
columns of the previously determined reference structure matrix are 
plotted. The direction numbers of the new long reference vectors are 
estimated. The direction cosines are determined and the fourth iterative 
stage reference structure matrix is computed from F. The fourth 
iterative stage is the final iterative stage for the illustrative 
example. The reference structure matrix and the reference vector 
intercorrelation matrix as determined by the fourth iterative are 
presented in Tables 13 and 14 respectively. 



Table 13 



Reference Structure Matrix As Determined 
By The Fourth Iterative Stage* ~ Matrix F 





A 4 

! 


B 4 


C 4 


1 


-.006 


.965 


.001- 


2 


.032 


.009 


.938 


3 


.964 


-.010 


.016 


4 


.025 


.462 


.723 


5 


.854 


.370 


-.002 


6 


.810 


-.009 


.387 


7 


-.003 


.707 


.510 


8 


.395 


.840 


.011 


9 


.468 


-.031 


.756 


10 


.015 


.602 


.621 


11 


.649 


. 654 


.002 


12 


.649 


-.027 


.605 


13 


-.003 


.951 


.032 


14 


.090 


-.040 


.919 


15 


.943 


.012 


-.026 


16 


.721 


.263 


.356 


17 


.294 


.629 


. 456 


18 


.016 


.943 


-.035 


19 


-.014 


.058 


.905 


20 


.955 


-.013 


.026 



*Thurstone, 1947, p. 213 



Table 14 

Intercorrelations of Unit Length Reference Vectors 
As Determined In The Fourth Iterative Stage* - Matrix Y 





A 4 


B 4 


C 4 




1.000 


-0.066 


-0.188 


B 4 


-0.066 


1.000 


-0.226 


C 4 


-0.188 


-0.226 


1.000 



*Thurstone, 1947, p. 213 
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Additional Algebraic Aspects of Thurstone* s Iterative Methodology: 

After the first iterative stage Thurstone defined, algebraically, a 
second type of transformation matrix H , This transformation matrix, 

H f - was used in conjunction with F m , the previously computed reference 



re- 



structure. matrix, to compute the u- th reference structure matrix, F„ . 

. w 

The matrix equations for computing F , disregarding H . , are reported by 

Equation 4. 

F = FV S D" 1 
u m mu u 

F u = FL Iu 



F = FV 
u u 



[4] 



Thurstone defined L in terms of 7 and S , equation 1. The matrix 

u m mu 

L is geometrically meaningful. However, H is defined in terms of S ' 

a inn uili 

-i 

and D and does not appear to be geometrically meaningful. 

H — c ! D~ ^ 

mu mu u 

The matrix F is defined as the product of the previous reference 

structure matrix, F , post-multiplied by . 

F - F H 
u m mu 



[5] 



F u = 



mu 



[ 6 ] 

[7] 



The matrix H therefore transforms F to F . The elements of R m . 

mu mu. mu 

are neither direction cosines nor direction numbers. The entries of 

TTili 

are expressed within the metric of the m ~ th stage unit reference vectors 

2 

while the entries of D u are expressed within the metric of the original 
fixed orthogonal frame. Any transformation matrix may be expressed as 
a product of all previous ^-matrices (Thurstone, 1947, p. 206). 



V u = (H 01 )(H 12 )(H 23 ) ( H s4 )...(H m ) 



[ 8 ] 
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For the illustrative example the //-matrices associated with third 
and fourth iterative stages are reported in Tables 15 and 16. For the 



30 



second iterative stage is identical to V ^ 

Table 15 

//-Matrix Computed For The Third 
Iterative Stage* - Matrix 11 

a 6 



0.945 


0.104 


0.229 


0.000 


1.042 


0.110 


-0.680 


-0.208 


0.917 



*Thurstone, 1947, p. 198 
Table 16 



//-Matrix Computed For The Fourth 
Iterative Stage* - Matrix //_ 

t54 



1.016 


.050 


.000 


.091 


.999 


.000 


.000 


-.050 


1.000 



*Thurstone, 1947, p. 209 

The function of may be thought of as providing an alternative 
method of computing F directly from F as opposed to computing F from 

W Tfl H 

F through the use of 7 . 

Summary of Section 1: 

Thurs tone's (1947) method of determining oblique transformations has 
been presented within the context of two dimensional sections. Through 
the use of an illustrative problem his terminology and matrices were 
discussed. Although certain aspects of Thurstone's approach were modified 
it may be assumed that the discussion presented in this section was taken 
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from Thurstone (1947, pp. 194-224) and typifies his approach to oblique 
transformation solutions. 

There are numerous objections to Thurstone' s procedures, all of 
which are either directly or indirectly associated with the matrix 5 . 

The matrix S is the subjective matrix of direction numbers defining 
the termini of the u- th iterative stage long reference vectors with 
respect to the unit length reference vectors of the ?.^-th iterative stage. 
It would be extremely difficult for two factor analysts working 
independently on the same factor matrix to determine identical oblique 
solutions when using Thurstone* s technique. It would be a most arduous 
task for a beginning student to apply Thurstone f s methodology successfully 
to a set of data in which (r > 3). Thurstone 's method becomes prohibitive 
timewise as the number of factors increases inasmuch as r(r - l)/2 plots 
are required at any one iterative stage to determine S . It is 
conceivable that a bad estimate of 5 might be obtained at some early 
iterative stage and not be recognized as such for several iterative 
stages. Finally, Thurstone 's method is just too time consuming and 
unreliable for all except the most experienced factor analyst. 

Section II 

A Basic Theoretical Rationale For The Simplified 
Obliquimax As Provided By The General Obliquimax 

In this section the general obliquimax is briefly discussed x^ith 

respect to the matrix equations defining an oblique solution. Special 

emphasis is placed on the matrix of direction numbers, L as discussed 

U 

in Section I, and the solution matrices expressed within the metric of 
the original factor solution. This section discusses oblique solutions . 
within the framework of variance modification and allocation. The 
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function of this section is to provide a basic theoretical rationale for 
the development of the simplified obliquimax in Section III. 

Common Variance Modification and Allocation: 

In determining an oblique transformation, a majority of the reference 
vectors, if not all of them, will covary having non-zero perpendicular 
projections upon each other as opposed to an orthogonal transformation in 
which none of the factor axes covary. Within the obliquimax framework 
the total common variance associated with an initial factor so _ution is . 
defined as the total sum of the squared projections of the variable 
vectors onto the initial factor axes. The common variance associated 
with any one of the initial common factors is defined as the sum of the 
squared perpendicular projections of the variable vectors onto that 
factor axis. An orthogonal transformation will not change the total 
common variance but it will in general alter the sum of the squared 
perpendicular projections associated with each common factor. 

When working within an oblique framework the process of column 
normalizing the direction numbers to form the matrix of direction cosines 
is actually a process of converting the metric from that of the initial 
fixed frame to the metric of the unit length reference vectors of the 
particular iterative stage associated with the direction cosines. Thus, 
for each iterative stage of an oblique solution the metric is changed' 
and is not comparable to the metric of the initial fixed frame. Because 
of this metric -variation the sum of the squared perpendicular projections 
of the variable vectors onto the reference vectors are neither comparable 
between iterative stages nor comparable with the sum of the squared 
projections onto the factor axes of the initial fixed frame. 
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The metric of the initial fixed frame may be retained in an oblique 
solution ii the matrix of direction numbers, previously defined by 
Equation 1, is used in place of the matrix of direction cosines, 
previously defined by Equation 3. The entries of the resulting 
"structure"* matrix would represent the perpendicular projections of the 
variable vectors onto the reference vectors, but the entries would be 
expressed within the metric of the initial fixed frame. Within 
Thurstone’s terminology such a matrix would be referred to as a long 
reference vector structure* matrix inasmuch as the entries are analogous 
to perpendicular projections of the variable vectors onto the long 
reference vectors as opposed to unit length reference vectors. When the 
metric is held constant in this fashion, it will be observed that the 
total sum of the squared projections of the variable vectors onto the 
long reference vectors is less than the total sum of the squared 
projections of the variable vectors onto the initial fixed axes. In the 
next subsection it will be demonstrated that the difference between the 
sum of the squared projections of the two matrices is accounted for by 
the perpendicular projections of the primary vectors onto each other, 
the covarying of the primary vectors. 

Within this framework the role of the direction numbers of the long 
reference vectors with respect to the initial fixed axes, I > , is one of 
defining the long reference vectors in such a manner that the sum of the 
squared perpendicular projections of the variable vectors onto these long 

* In a latter portion of the next subsection it will be demonstrated 
that the long reference vector structure matrix is not the only meaningful 
name that might be given to the matrix being discussed. 
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reference vectors is less than the sum of the squares of their projections 
onto the initial factor axes. There should be some systematic relationship 
between the matrices of direction numbers determined in successive 
iterative stages such that the sum of the squared perpendicul ir projections 
of the variable vectors onto the long reference vectors becomes successively 
smaller at each stage. It will be demonstrated that such a framework will 
provide numerous conceptual and algebraic advantages. 

Theoretical Equations for the General Oh liquimax Tra?is formation: 

The objective of this subsection is to provide a basis for the 
computational aspects of the simplified obliquimax through a brief 
presentation of the equations and logic for the general obliquimax (Hofmann, 1971). 

Assume some initial factor loading matrix F, defining the perpendicular 
projections of n variable vectors onto r mutually orthogonal factor axes as 
determined by the initial factoring method. The problem is to select by u, 
successive approximations the unit reference vectors A , B. and C such 

ui (t ui 

that the number of variable vectors with vanishing projections onto these 
unit reference vectors is a maximum. (Where deemed necessary for 
illustrative purposes r will be assumed to be three.) 

In the obliquimax transformation all matrices of direction numbers 

* 

0 

are defined as the product of some positive definite diagonal matrix and 
some orthonormal transformation matrix, T. A matrix of direction numbers 
developed in this manner, and the ensuing matrix of direction cosines, will 
always be non-singular and generally non-orthogonal. This somewhat unusual 
approach was first suggested by Harris and Kaiser (1964) in their classic 
paper on determining oblique transformation solutions through the use of 
orthogonal transformation matrices. Although the discussion in this paper 
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will only encompass their case I and case II solutions, the general 
obliquimax encompasses all three of the cases discussed by Harris and 
Kaiser. 

t 

o 

Let the ii-th element of the positive definite diagonal matrix p w , 
represent the sum of the squared projections of the r, variable vectors 
onto the i- th factor axis associated with F. 

D % = diagonal [F’F] [9] 

All matrices of direction numbers within the general obliquimax 
framework are defined specifically as the product of some exponential 
power, p, of D and an orthonorraal matrix T. Equation 10 represents the 
matrix of direction numbers, L , for defining the termini of the w-th 
iterative stage long reference vectors with respect to the initial fixed 
f r ame . 

L u « TT^T [10] 

Along with the algebraic advantages of such a definition of L , there 
are several conceptual and interpretative advantages. In defining the 
positive definite diagonal matrix as some exponential function of the 
column sums of squares of F , the matrix of direction numbers is implicitly 
some function of the common variance and hence a function of the 
perpendicular projections of the variable vectors associated with F. 

The matrix F may be rewritten as: 

F = ED. [11] 

The matrix E is the column normalized form of F and the elements of D 

2 

are the square roots of the elements of D . The matrix F may be post 

£ 

multiplied by L.. to form F which would be the matrix whose elements 
represent the perpendicular projections of the variable vectors onto the 
long reference vectors determined by the w-th iterative stage. 

F u = FL u = fd ~ Pt 

o 

ERJC . 
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There are several important aspects of Equation 12 that should be 
noted. Any matrix of direction numbers within the general obliqulraix. 
framework is simply a row rescaling of T by some exponential power, r, 
of the matrix D. Equation 12 may be rewritten as Equation 13. 

F* = EDD~ V T [13] 

1st 

The discussion of common variance associated with ?*, will be based unon 

It 

the sum of the squared perpendicular projections onto the r long reference 
vectors associated with F . Inasmuch as the columns of E are normalized, 
the matrix product (DD &T) defines the allocation of the perpendicular 
projections of the variable vectors onto the long reference vectors. That 
is to say, the common variance, within the metric of the initial fixed 
frame, associated with F can be discussed with respect to the matrix 

1st 

(1 DD~ P T ) . 



The matrix T is an orthonormal matrix o.f direction cosines. Therefore, 
the sum of the squares of any column or row of T is unity (T f T = TT r ~ I). 
The square of any cosine is a proportion. The square of any diagonal 
element of (DD~P) is a portion of the total common variance associated 

a 

with F. The variance aspects of F' can be discussed through a consideration 
of the squared elements of (DD ^T)> The variance discussed within this 
specific framework would be with respect to the sums of the squared 
perpendicular projections of the n variable vectors onto the r long 
reference vectors. A symbolic representation of the squared elements of 
(DD ^ T ) is reported in Table 17. 

The variance contributed by fixed axis A to the sum of the squared 

perpendicular projections of the n variable vectors onto all r long 

2 2r> 

reference vectors is the variance contributed by fixed axis 



, 

22 



is (d a ney /d 2 ^); 



the variance contributed by fixed axis C Q 
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Table 17 

Symbolic Representation of the Squared 
Elements of (DD"^T) (assume v = 3) 




(Other than noting that the exponent p will be chosen such that the total 

% 

column sums of squares of F is less than the total associated with F 

discussion of p will be deferred at this point.) The proportion of that 

variance contributed by A Qf (d^/d^), that will be allocated specifically 

to the sum of the squared projections of the variable vectors onto long 

2 2 2 

reference vector: A* is cos 6--; BJ, is cos 6„ 0 ; Cl is cos A 

w 11 ** 22 u 13 

similar interpretation may be made for the other (r - 1) fixed axes with 

respect to the sum of the squared projections of the n variable vectors 

onto the r long reference vectors. 

A portion of the total variance associated with a fixed axis is 

allocated to the sum of the squared perpendicular ‘projections of the 

variable vectors onto the long reference vectors. The discussion here 

concerns the proportionate distribution of this variance by a fixed axis 

to each of the long reference vectors. The portion of variance associated 

with the j-th fixed axis that is allocated to the sum of the squared 

perpendicular projections of the n variable vectors onto the r long 

2 2v 

reference vectors is always (djj/djj ) . Of that particular portion of 

<7 d vd 

variance, the proportionate allocation to the i-th long reference vector 
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is Q 08 (i « 1, 2, . . .r) . Therefore, at any particular iterative 

stage the only value that will be altered in Table 17 is the exponent p . 
(In the definitive paper on the general obliquimax, Table 17 is further 
modified to provide additional information about the reference vectors.) 

To further clarify the role of variance modification in the general 
obliquimax, it is necessary to discuss the variance covariance matrices 
associated with both the reference vectors and the primary vectors. Let 
the symmetric matrix i?*, of order n by w and singular of rank **, be the 
major product of F . 

R* = FF* [13] 

The matrix R* is a reproduced matrix whose off-diagonal elements 
approximate either the observed correlations or covariances between the 
n variables. 



For any oblique solution Equation 14 must hold. 

R* = F (Y )~^F' 

n u [JL u J u 



The matrices F u and Y u were previously defined by Equations j* and 5. 



A 



[14] 



Let the matrix Z u represent the matrix of intercorrelations associated 

with the r unit length primary vectors. Then by definition, Z u is the 

2 

normalized inverse of Y u (Thurstone, 1947, p. 215). Let the matrix D ^ be 



,-l 



defined as the diagonal of Y , therefore pre- and post-multiplication of 
-1 - 1 

Y u by Dg will result in the matrix of intercorrelations for the primaries. 



= diagonal [J“^] 

Z = D~ 1 (Y~ 1 )D~ 1 

u 3 K u J 3 



[15] 

[16] 



Through substitution from Equations 2, 3, 5 and 10, the matrix Y may be 
defined algebraically within the obliquimax framework. 

Y u = D^ 1 T , D’ 8p TD^ 1 [17] 



o 
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Following from Equation 17 is the algebraic definition of within the 
obliquimax framework. 

Z u = DZ 1 (D,.T , D 2p TD u )D’l 1 = CrV,) (T'D 2p T) (DDZ*) [18] 

w 6 u u 3 o u « 3 



- 1 . 



In Equation 18, the matrix product of (£7 is just the diagonal mitrix 



n __ x 

that normalizes the columns of T ) . The matrix (L/L.Y Is the variance 
covariance matrix, Z , associated with the long primary vectors. Thus, Z 



-7 



M 5 



is just the inverse of Y u > 

Y* = T'D~ 2p T 

u 



Z* = T’D Sp T = Cl *; -2 



[19] 

W - - “ ’-IT ^0] 

It is immediately apparent from Equations 19 and 20 that the only 

algebriac difference between the primary variance covariance matrix and 
the reference variance covariance matrix is the sign of the exponent p. 
Furthermore, the only algebraic difference between the matrix of direction 
numbers, (D P T ) , for the primary structure matrix, and the matrix of 
direction numbers, ( D P T) , for the reference structure matrix, is the sign 
of the exponent p. Thus, the matrix of direction numbers for the long 
primary structure matrix is just the transpose of the inverse of the 
matrix of direction numbers for the long reference structure matrix. 
Equations 19 arid 20 incorporate the metric of the initial fixed axes. 
Utilizing Equations 19 and 12 it is possible to redefine R* using 



matrices expressed in the metric of the initial fixed axes. 
R* = F* (Y*)~^F* r 

By substituting % u for Equation 22 results. 

= F X F u' 



[ 21 ] 



[ 22 ] 



The trace of R* represents the total common variance and it is equal to 
the total column sums of squares for the matrix F> Equation 22 is presented 
to demonstrate that the total variance associated with F may be thought of as 
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being distributed between F* and Z*. In explaining Table 17 the variance 

j 

associated with F?. was discussed with respect to the sum of the squared 

perpendicular projections of the n variable vectors onto the v long 

reference vectors. It was pointed out that the variance associated with 

F* would always be less than the variance associated with f. It may 

# 

therefore be concluded that the common variance not allocated to F u must 
be allocated to Z*. The variance allocated to Z^ accounts for the 
perpendicular projections of the 2? long primary vectors onto each other. 

The variance associated with Z^ can be discussed with respect to (jPt) , 
the matrix of direction numbers for computing the primary structure 
matrix. Although interest here does not center on the primary structure 
matrix, it may be inferred from Equation 22 that the matrix of direction 
numbers associated with such a matrix is instrumental in explaining 
variance allocation. The variance aspects of Z u can be discussed through 
a consideration of the squared elements of (Lp T ) . A symbolic representation 
of these squared elements is reported in Table 18. 

Table 18 



Symbolic Representation of the Squared 
Elements of (JpT) (assume r = 3) 



A 9 - 

u 


-Primary 


B^~ Primary 


c u- 


-Primary 


d 2 P 

11 


2 * 
cos 3 


nOV 2^ 

d il ° 0S >12 


d 2p 

11 




£v 

22 


2 o 

C0S hi 


Jp 2 - 

d 22 003 0 Si 


Jp 

U ^ 

i-j »'• 


°°*\z 


d 2 £ 

33 


2 n 

cos $• 

31 


eft cos 2 3 

33 32 


eft 

33 


2. 

COS (J 

•7 7 
0 



The variance contributed by fixed axis A Q to the sum of the squared 
perpendicular projections of the long primaries onto each other is (cLO; 
the variance contributed by fixed axis B Q is (< $22 ) » the variance 
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contributed by fixed axis C is (d^) . The proportion of that variance 

^ 66 



contributed by A Qi u i r/), that will be allocated to the sum of the squared 



perpendicular projections of the ( 2 ? - 1) long primary vectors onto long 



interpretation may be made for the other (r - 1) fixed axes with respect 
to the sum of the squared perpendicular projections of the r long primary 
vectors onto each other. * . 

A portion of the total variance associated xd.th a fixed axis is 
allocated to the sum of the squared perpendicular projections of the long 
primary vectors onto each other. The discussion with respect to Table 18 
concerns the proportionate distribution of this variance by a fixed axis 
to each of the long primary vectors. The portion of variance associated 
with the j-th fixed axis that is allocated to the sum of the squared 
perpendicular projectims of the r long primary vectors onto each other 
is always (d^£) . Of that particular portion of variance, the proportionate 

allocation to the i-t\\ long primary vector is cos (£ = 1, 2, ... 2 ?). 

J'Z' 

Therefore, at any particular iterative stage, the only value that will be 

altered in Table 18 is the exponent p. 

Comparing Tables 17 and 18, several generalizations may be made. 

2 

The variance associated with the 'k-th fixed axis, d.^, may be divided 

2 2v 

into two multiplicative portions (d^/d^s) and (d^) • One portion of* the 
2 2 v 

variance, (d. »/d..). is associated with the sum of the squared 
perpendicular projections of the n variable vectors onto the r long 
reference vectors. The second portion of the variance, (d|?) , is 
associated with the sum of the squared perpendicular projections of the 
r long primary vectors onto each other. The proportionate contribution 



2 2 

primary vector: A ^ is cos 3 B ^ is cos is co 








made to the sum of the perpendicular projections of the tl 



ERIC 



*The element £• cos 3„* j also represents* the perpendicular projection of 
the long reference vector 3 onto the i-th fixed axis. 
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variable vectors onto the J-th long reference vector is cos 3 • • . The 

I'd 

proportionate contribution of (dfo) made to the sum of the squared 

perpendicular projections of the (r - 1) long primary vectors onto the 

2 

j-th long primary vector is cos 3* *. Equation 23 defines the division 
of variance for the £-th fixed axis. 

4 =( 4 / 4 >< 4 >- (233 

Several observations with respect to p can be made from Equation 23. 
If (p = 0), there will be no allocation of variance to the projections of 
the primaries onto each other and the associated solution will be an 
orthogonal transformation solution. If (p = 1) , all of the variance will 
be allocated to the projections of the primaries onto each other and the 
resulting oblique solution would be analogous to the Harris and Kaiser 
(1964) independent cluster solution. As the value of p progresses from 
zero to unity, the amount of variance allocated to the projections of the 
primaries onto each other becomes progressively larger. Thus, as p 
progresses from zero to unity, one should expect the primary vectors to 
become progressively more correlated. One might be inclined to limit the 
values of p to values within the interval bounded by zero and unity, but 
there is at this time no rationale for such a limitation. (For some sets 
of empirical data the general obliquimax has iterated to a value of p 
that is slightly larger than unity.) It is prudent to realize that for a 
value of p considerably larger than unity, but not necessarily larger 

-p 

than three, some of the diagonal entries .of D c might approach zero and 

-p 

the matrix D will not for all practical purposes be a positive definite 
matrix. (Theoretically, D ^ will always be positive definite however 
computationally some near zero values will function as zeros.) 
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Within the framework defined, the only algebraic value that changes 
from one iterative stage to the next in defining ? u is the exponent p. 

The discussion, however, has centered around the long reference vectors 
and long primary vectors. Traditionally an oblique solution is interpreted 
within the framework of either unit length reference vectors or unit 
length primary vectors. In order to establish the oblique solution within 
a traditional framework it will be necessary to compute the diagonal 
matrices that will rescale F u , Y u and Z u to the metric of unit length 
vectors. 

From Equations 2, 3 and 4 it may be inferred that the diagonal matrix 

for rescaling F*. to the unit length reference structure matrix is 

u 

determined from the diagonal of (L^L^) . In the previously mentioned 

_2 

equations this diagonal matrix was referred to algebraically as D . 

Within the obliquimax framework this diagonal matrix will be referred to 

-1 2 , 

as D u The matrix D ^ is defined by Equation 24. 



D 2 u1 = diagonal [T'D~ Sp T] 

Equation 4 may be rewritten as Equation 25 to define the unit length 

reference structure matrix within the algebraic framework of the 

obliquimax. . 

F = FD~ P TD~l 
u ul 



[24] 



[25] 



In Equation 20 the variance covariance matrix of the long primary 
vectors, Z*, was defined as the inverse of the long reference vector 
variance covariance matrix. It was inferred from Equations 19 and 20 
that the matrix of direction numbers for transforming F to the long 
primary structure matrix was just the transpose of the inverse of the 
matrix of direction numbers necessary to transform F to the long 
•reference structure matrix. However, the primary structure matrix is not 



o 
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generally used for interpretative purposes. It is the primary pattern 
matrix that is interpreted within Holzinger’s framework. The matrix 
necessary to transform the initial matrix F to the primary pattern 
matrix, denoted as W f is by definition the transpose of the inverse of 

(A 

the primary structure transformation matrix. The primary structure 

transformation matrix is just the column normalized form of the 
-1 

direction numbers associated with the long primary vectors. Let D . 
represent the diagonal matrix that normalizes the columns of (D^T) . 

Then from Equations 24 and 25 Equations 27 and 28 follow as: 

D u2 = dia 9 onal \.T'D 2p T ] ; 

z nj = d~It’d 2 Ptd~1. 

u u2 u2 

Let the matrix (JPTD represent the transformation matrix for the 
primary structure matrix. The transpose of the inverse of this matrix is 

The matrix necessary to transform F to the primary pattern matrix, 

— n 1 

W ui as determined by the w-th iterative stage is (Z? ^TD ) . 

w.. = fd~ p td 2 0 

u u2 



[27] 

[28] 



[30] 



However, if reference is made to the parallel projections of the n 
variable vectors onto the r long primary vectors, then the entries being 
referred to are parallel projections within the metric of the initial 
fixed frame. That is to say, W* is just W without the column rescaling. 

[31] 

Comparing Equation 31 with Equation 12 it becomes immediately evident that 

[32] 

The Holsinger parallel projections of the n variable vectors onto 
the r long primary vectors cere identical to the Thurstone perpendicular 
projections of the n variable vectors onto the r long reference vectors . 



J7* = FD~ P T 



w u = F u = FD ' ?T - 
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That is to say 3 when the metric of the initial factor solution is retained 
in place of either unit length reference vectors or unit length primary 
vectors the Holzinger pattern matrix and the Thurstone structure matrix 
are identical. (In earlier papers on the obliquimax W u and F u were 
referred to as one matrix, the "basic matrix.") 

The above paragraph follows logically from an earlier paper presented 
by Harris and Knoell (1948) in which they derive a diagonal matrix that 
will rescale the columns of a Holzinger pattern matrix to those of a 
Thurstone structure matrix. The inverse of this matrix will rescale the 
columns of a Thurstone structure matrix to those of a Holzinger pattern 
matrix. Their discussion centered around the geometry of the two 
solutions and was considerably less complex, algebraically, than this 
section. 

In their discussion Harris and Knoell (1948) demonstrate that the 
Thurstone structure values represent the bases and the Holzinger pattern 
values represent the hypotenuses of similar right-angled triangles. The 
elements of the diagonal rescaling matrix that they use for their 
conversion represent the correlations between the reference vectors and 
their associated primary vectors. They note that a primary vector is 
defined by the intersection of (r - 1) hyperplanes and it is uncorrelat id 
with the normals to these hyperplanes. Therefore, each primary vector is 
by definition orthogonal to all but one reference vector. Each reference 
vector is correlated with just one primary vector. 

The correlation between a unit length reference vector and a unit 
length primary vector is defined as the scalar product of the paired 
direction cosines of the vectors. Equation 33 defines the diagonal 
matrix £> whose ii-th element represents the correlation between the i-tb. 

o 
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[ 33 ] 



reference vector and its associated i-th. primary vector. 



D * - (D-^T'D-P) = D- u ]D- u l 



Swrnary of Section II 

The basic framework for the general obliquimax has been discussed. 
The algebraic definitions of certain oblique solutions were discussed 
within the framework of common variance and perpendicular projections. 
The "solution equation^' of the general obliquimax may now be summarized. 



— P — 1 

Reference Structure = FD 

Primary Pattern = FD~^TD U 2 

Primary Intercorrelations - dJ^T ' D^TD^ 

- 1 - 2‘0 ~2 

Reference Intercorrelations = D ul T'D -ro* 

Inter correlations Between Primaries and Reference Vectors 



n -l -2 
D ul D u2 



The discussion presented in this section is by no means complete. 
There has been a minimum of discussion concerning the exponent p and 
there has been no discussion with respect to the specific computation of 
T . In the next section the computation of T will be briefly discussed 
and the use of the exponent p will be somewhat clarified. This section 
has provided a basic rationale for the equations to be used in the 
simplified obliquimax. 



Section III 

The Simplified Obliquimax 

Having discussed certain interpretative properties of the direction 
numbers and defined in part the analytic computations of the direction 
numbers it is now possible to re-establish the problem of computing an 
oblique transformation within an obliquimax framework. The iterative 
procedure developed in this section, referred to as the simplified 
‘obliquimax, is developed for pedagogical purposes. Certain matrices are 
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modified in the iterative process to provide computational simplicity. 

The simplified obliquimax is a reliable semi-subjective transformation 
procedure that will allow beginning students to successfully compute an 
oblique transformation solution. It is also hoped that the simplified 
obliquimax will clarify certain aspects of oblique transformations in 
general. 

Assume some factor loading matrix, F, defining the perpendicular 
projections of n variable vectors onto v mutually orthogonal factor axes. 

The p axes are arbitrarily orthogonal axes as determined by the initial 
factoring method. For illustrative and comparative purposes the 
simplified obliquimax will be discussed within the framework of 
Thurstone’s (1947, pp. 140-144) classic box problem. The centroid 
solution for the box problem is reported in Table 1. For this particular 
set of data (n = 20) and (p « 3) . 

The P axes may be denoted as A 0 , B 0 and C Q . These arbitrary axes 
are regarded as fixed in position. The problem is to select by u 
successive approximations the unit reference vectors, A u , B u and Cu> 
such that the number of variable vectors with vanishing projections onto 
these unit reference vectors is a maximum. 

With the exception of the exponent p and the final diagonal rescaling 
matrices all matrices necessary for computing either a reference structure 
or a primary pattern matrix can be determined prior to the computation 
• of p. Thus, the problem at hand is one of determining a value for the * 

exponent p such that the associated reference structure matrix will have 
a maximum number of vanishing projections. 

In starting the iterations for the simplified obliquimax the matrix 

o 

D2 is determined from F through the use of Equation 9. Presented in 

o 

Table 19 is whose ii - th element represents the total sum of the 

o 
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squared projections of the n variable vectors onto the i- th fixed axis. 

= diagonal [F'F] [ 9 ] 

Table 19 

The Positive Definite Diagonal Matrix Whose ii-th 
Element Represents the Variance Associated with 
the i- th Fixed Axis ~ Matrix D * 





B o 




12.561 


.000 


.000 


.000 


3.920 


.000 


.000 


.000 


3.154 



The orthonorraal transformation matrix T used in the simplified 
obliquimax is computed from F using Kaiser’s (1958) normal varimax 
transformation procedure. A complete rationale for the use of the normal 
varimax transformation is beyond the scope of P'is paper. A very basic, 
though not particularly compelling, rationale for its use is that a 
primary vector may be defined in a special sense as a linear least squares 
approximation to a group of variable vectors without the restriction of 
mutual orthogonality and the normal varimax axes are with the restriction 
of mutual orthogonality in a special sense a linear least squares 
approximation to a group of variable vectors. It is assumed, and as will 
be seen, this is a critical assumption, that if the restriction of 
orthogonality were placed on the primary or reference vectors they would 
occupy the precise positions of the varimax axes. 

Initially the matrix F must be transformed orthogonally to determine 
T. The matrix FT will be utilized in the simplified obliquimax as opposed 
to F. The matrix FT may be thought of as a preparatory matrix for the 
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iterative stages (see Table 20 for FT) . The varimax axes will be assumed 

as the initial V unit length reference vectors and will be referred to as 

A v , B v and C v . (The subscript V referring to orthogonal varimax axes.) 

# 

As in Thurstone’s procedure the v unit length reference vectors are 

initially assumed to be mutually orthogonal. Contrary to Thurstone's 

0 

procedure the initial unit reference vectors of the simplified 
obliquiraax are not collinear with the r initial factor axes. 

Table 20 

Initial Factor Matrix Orthogonally 
Transformed By Normal Varimax - Matrix FT 





Ay 


B v 




1 


0.052 


-0.990 


-0.109 


2 


0.146 


-0.142 


-0.973 


3 


0.984 


-0.050 


-0.103 


4 


0.135 


-0.576 


-0.803 


5 


0.887 


-0.431 


-0.117 


6 


0.872 


-0.094 


-0.472 


7 


0.091 


-0.797 


-0.605 


8 


0.443 


-0.887 


-0.141 


9 


0.566 


-0.102 


-0.820 


10 


0.118 


-0.705 


-0.710 


11 


0.692 


-0.711 


-0.134 


12 


0.733 


-0.095 


-0.681 


13 


0.046 


-0.980 


-0.139 


14 


0.200 


-0.092 


-0.954 


15 


0.958 


-0.065 


-0.060 


16 


0.790 


-0.362 


-0.462 


17 


0.384 


-0.726 


-0.568 


18 


0.057 


-0.963 


-0.071 


19 


0.097 


-0.185 


-0.940 


20 


0.975 


-0.048 


-0.112 



Although graphical plots are not necessary in computing the simplified 
obliquimax they are presented initially for comparative and illustrative 
purposes. Figures 7, 8, and 9 represent the planar plots of the variable 
points with respect to the fixed axes. Also included in these figures are 
the projections of the unit reference vectors Ay, By and Cy. 
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Figure 7 

Planar Plots of Variable Vector Termini with Respect 
to Fixed Axes A Q and B Q9 Projected Onto Plane A Q B Q 




Figure 8 

Planar Plots of Variable Vector Termini with Respect 
to Fixed Axes A Q and C Q , Projected Onto Plane A Q C 0 
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Figure 9 

Planar Plots of Variable Vector Termini with Respect 
to Fixed Axes B Q and C Q , Projected Onto Plane B Q C Q 




The normal varimax transformation matrix, T> is a matrix of direction 
cosines defining the unit reference vectors A v> By and C v with respect to 
the fixed axes A q9 B q , C Q . This matrix will be utilized in all iterative 
stages and is presented in Table 21. 

Table 21 

The Orthonormal Transformation Matrix Computed From F By 
The Norma"?. Varimax Transformation Procedure - Matrix T 





1 A v 


V 


c v 


A o 


.587 

* 


-.572 


-.574 




.564 


.797 


-.217 


V 


.581 


-.196 


.790 
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In the previous section it was established that a good estimate of 
p will usually be some number bounded by zero and unity. Select some 
small number less than unity but greater than zero and call that number 
x. The initial estimate of the exponent p will be x* For the 
illustrative example x is chosen to be .20. (This choice for the 
illustrative example is not arbitrary. The rationale for choosing .20 
instead of .05 or .50 will become apparent. In practice a choice of .10 
would be a very safe estimate.) 

Using x , T and D a matrix of direction numbers is defined as T r D T. 

This matrix of direction numbers defines the termini of a new set of long 
reference vectors with respect to the axes of the orthogonally 
transformed matrix FT or the varimax axes A v , B v and C v . Because the 
matrix T is orthonormal post-multiplying FT by T , D'~ X T will result in an 
Equation, 34, analogous to Equation 12. 

[ 12 ] 

[34] 

Effectively, the post-multiplication of (FT) by (T f B~ X T) is the same as 
post-multiplying ( F ) by (D~^T) if p = x» The matrix (T’D ~T) corresponds 
to the subjective Thurstonian matrix S at this point of the paper, 

JL u 

however it will be referred to simply as S and for subsequent iterative 
stages it will be somewhat different in function than Thurstone’s S mu matrix. 
S = T’D~ X T * * ' [35] 



u 

,* 



F’u = FD~Vt 



F* = (FT)(T'D~ X T) = FD~ X T 

1 



In Figures 1C), 11 and 12 the variable points have been plotted with respect 

i 

to the columns of FT (see Tabl*?. 20). The termini of the long reference 
vectors Al , , B’ and C 1 have been plotted with respect to A v , B v and <7 , 
the varimax axes. The coordinates for the termini of the new long reference 

4* 0 

vectors with respecv. 1 to varimax axes A v , B v and C v (initial r unit length 
reference vectors) are given by the columns of the matrix S in Table 22. 



o 

ERIC 



53 



* 



t 



» 



53 



Table 22 

Symmetric Matrix of Direction Numbers of First 
Iterative Stage Reference Vectors with Respect 
To Previous Reference Vectors* - Matrix S 





A f 

A u 


B' 

u 


rt 

U U 


A' , 
u-1 


.846 


.030 


.041 


B 1 

U-1 


. .030 


.842 


-.034 


n t 

"u-1 


.041 


-.034 


.853 



^Subscripts will be referred to in a 
later portion of this paper. 

Figure 10 

Planar Plots of Variable Vector Termini with Respect 
To Varimax Axes A v and B v , Projected Onto Plane A^B^ 
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Figure 11 

Planar Plots of Variable Vector Termini with Respect 
To Varimax Axes Ay and Cy, Projected Onto Plane Ay 7y 




Planar Plots of Variable Vector Termini with Respect 
To Varimax Axes By and Cy, Projected Onto Plane ByC v 
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Equation 34 may be rewritten as Equation 36 to specifically define 
the long reference structure matrix with respect to the first iterative 
stage. * 

F* = .FTS = FD~ X T ' [36] 
The first iterative stage long reference structure matrix is reported in 
Table 23. 



Table 23 

Long Reference Structure Matrix as Determined 
By First Iterative Stage - Matrix Fj 





A r 

R 1 


B 'i 


c i 


1 


0.010 


-0.828 


-0.057 


2 


0.079 


-0.082 


-0.819 


3 


0.826 


-0.009 


-0.046 


4 


0.064 


-0.454 


-0.659 


5 


0.733 


-0.332 


-0.048 


6 


0.715 


-0.037 


-0.364 


7 


0.028 


-0.647 


-0.485 


8 


0.343 


-0.728 


-0.072 


9 


0.442 


-0.041 


-0.673 


10 


0.050 


-0.566 


-0.577 


11 


0.559 


-0.573 


-0.061 


12 


0.589 


-0.035 


-0.547 


13 


0.004 


-0.819 


-0.083 


14 


0.127 


-0.039 


-0.802 


15 


• 0.806 


-0.024 


-0.010 


16 


0.638 


-0.265 


-0.349 


17 


0.280 


-0.580 


-0.443 


18 


0.017 


-0.807 


-0.025 


19 


0.038 


-0.120 


-0.791 


20 


0.819 


-0.007 


-0.054 



Within Thurs tone's algebraic framework the initial matrix of 
direction cosines, V Q , is post-multiplied by the matrix S 0 j to form the 
matrix The matrix Lj is the matrix of direction numbers for the long 

reference vectors with respect to the fixed frame. Within the 
obliquimax framework V Q is just T and S Q j is S, therefore: 



*It is important to note here that unlike Thurstone (1947) we do not 
assume F to represent the first iterative stage of the oblique transformation 
0 •" solution, nor do we assume FT to represent the first iterative stage. 
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L 1 = V S - = T(.T'D~ X T) = ZT^. • [37] 

i o ol 

If a; is equal to p Equation 37 would be identical to Equation 10. 

l u = zrPr [io] 

The variance-covariance matrix of the first iterative stage reference 
% 

vectors, is computed using Equation 19 and may be inferred from 
Equation 5 and 37 as: 

Yj = T’D~ 2X T = S 2 = SS = (T'D~ X T) (T'D~ X T) = 5 '5. [38] 

Inasmuch as the matrix of direction numbers, 5, is symmetric it is 

* 2 

sufficient to simply multiply it by itself and note that Y = S . 

Table 24 



Variance Covariance Matrix of the First Iterative 
Stage Long Reference Vectors - Matrix Y* 





Al 

i 


B i 


f" 


A i 


0.713 


0.049 


0.069 


B 'l 


0.049 


0.711 


-0.057 


C 1 


0.069 


-0.057 


0.730 



The variance-covariance matrix, (see Table 25) , of the first iterative 

stage long primary vectors, Z^ , is just the inverse of Y^ , (see Table 24), 

-1 

which is the product of the inverse of the direction numbers, S » 

_ I 

pre-multiplied by its transpose, (5 )'• Inasmuch as S is a symmetric 

matrix Z* is just S 2 . 

Z* = T’D 2x T = S~ 2 = S~ 1 S~ 1 = (JT’D X T' (Y"D X T) = (Y*)" 2 [39] 

The problem of coalescing of the reference vectors has been 
eliminated by the definition of the obliquimax direction numbers. 

Therefore, unlike Thurstone’s solution it is not necessary to rescale 
Y 2 to form Y ' to check for coalescing of reference vectors. At this 
point of the discussion it will be assumed that a second iteration is 
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Table 25 

Variance Covariance Matrix of the First Iterative 
Stage Long Primary Vectors - Matrix Z* 





A 9 -primary 


B ^-primary 


C'-vrimary 
1 * 


A* -primary 


1.416 


-0.110 


-0.140 


B* -primary 


-0.110 


1.427 


0.123 


C ^-primary 


-0.142 


0.123 


1.395 



needed. A discussion of the evaluation of the long reference vector 
structure matrix will be deferred until after the discussion of the 
computation of subsequent iterative stages. (For the reader’s interest 
the first iterative stage unit length reference structure matrix and unit 
length reference vector intercorrelation matrix are reported in Tables 26 
and 27 respectively.) 




Table 26 



Reference Structure Matrix as Determined 
By The First Iterative Stage - Matrix F 







B i 


c i 


1 


0.012 


-0.982 


-0.067 


2 


0.093 


-0.097 


-0.958 


3 


0.975 


-0.011 


-0.054 


4 


0., 075 


-0.538 


-0.772 


5 


0.865 


-0.394 


-0.056 


6 


0.844 


-0.044 


-0,426 


7 


0.033 


-0.768 


-0.568 


8 


0.404 


-0.864 


-0.084 


9 


0-. 522 


-0.048 


-0.787 


10 


0.059 


-0.671 


-0.675 


11 


0.659 


-0.680 


-0.071 


12 


0.695 


-0.041 


-0.640 


13 


0.004 


-0.972 


-0.097 


14 


0.150 


-0.046 


-0.938 


15 


0.951 


-0.029 


-0.011 


16 


0.753 


-0.315 


-0.409 


17 


0.330 


-0.688 


-0.519 


18 


0.020 


-0.957 


-0.029 


19 


0.045 


-0.142 


-0.926 


20 


0.967 


-0.009 


-0.063 
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Table 27 

Intercorrelations of the Unit Length Reference Vectors 
As Determined By the First Iterative Stage - Matrix Y ^ 





A 1 




G i 




1.000 


0.069 


0.095 


B 1 


0.069 


1.000 


-0.079 


C 1 


0.095 


-0.079 


1.000 



A symmetric matrix of direction numbers, 5, was defined in the first 
iterative stage. This matrix of direction numbers will be a constant 
throughout all iterative stages. Specifically it describes the termini 
of the present iterative stage long reference vectors with respect to the 
previous stage long reference vectors. Here it must be emphasized that 
the obliquimax matrix of direction numbers S is distinctly different from 
Thurstone’s matrix of direction numbers S . For the first iterative 

ifll/i 

stage these two matrices have identical roles. On subsequent iterative 
stages the obliquimax matrix S defines the termini of the present iterative 
stage long reference vectors with respect to the previous iterative stage 
long reference vectors. On subsequent iterative stages the Thurstone 
matrix S mu defines the termini of the present iterative stage long 
reference vectors, u , with respect to the previous iterative stage unit 
length reference vectors , m. To compute the second iterative stage long 

JL 

reference structure matrix it is only necessary to post-multiply F „ by S . 

(See Table 28 for F*; for continuity F is reported in Table 29.) 

F* = FD~ 2X T = ( FD~ X T)(T'D~ X T ) = F*S [40] 

6 1 
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Table 28 

Long Reference Structure Matrix As Determined 
By Second Iterative Stage - Matrix F* 0 





A ' 

A 2 


B 2 


C 2 


1 


-0.019 


-0.695 


-0.020 


2 


0.031 


-0.038 


-0.692 


3 


0.697 


0.019 


-0.005 


4 


0.013 


-0.357 


-0.544 


5 


0.608 


-0.256 


0.000 


6 


0.589 


0.003 


-0.279 


7 


-0.015 


-0.527 


-0.390 


8 


0.265 


-0.600 


-0.022 


9 


0.345 


0.002 


-0.554 


10 


0.001 


-0.455 


-0.470 


11 


0.453 


-0.463 


-0.009 


12 


0.475 


0.007 


-0.441 


13 


-0.025 


-0.686 


-0.042 


14 


0.073 


-0.001 


-0.677 


15 


0.680 


0.004 


0.026 


16 


0.518 


-0.192 


-0.262 


17 


0.201 


-0.465 


-0.347 


18 


-0.011 


-0.678 


0.007 


19 


-0.004 


-0.073 


-0.669 


20 


0.690 


0.020 


-0.012 
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Table 29 

Reference Structure Matrix As Determined By 
Second Iterative Stage - Matrix F 0 



A 2 B 2 C 2 



1 


i— — 

; -0.026 


-0.972 


-0.027 


2 


| 0.042 


-0.054 


-0.941 


3 


i 0.964 


0.026 


-0.007 


4 


; 0.018 


-0.500 


-0.740 


5 


! 0.841 


-0.358 


0.001 


6 


1 0.814 


0.004 


-0.380 


7 


-0.021 


-0.738 


-0.530 


8 


0.367 


-0.840 


-0.030 


9 i 


| 0.477 


0.003 


- 0.753 


10 


0.002 


-0.636 


-0.639 


11 


1 0.626 


-0.648 


-0.013 


12 | 


0.657 


0.010 


-0.600 


13 


| -0.034 


-0.961 


-0.058 


14 


[ 0.101 


-0.002 


-0.921 


15 


1 0.941 


0.006 


0.035 


16 | 


0.716 


-0.269 


-0.357 


17 ! 


0.278 


-0.650 


-0.471 


18 


-0.015 


-0.948 


0.010 


19 


-0.005 


-0.102 


-0.910 


20 


0.955 


0.028 


-0.016 



Note that effectively (TT r = I) and that post-multiplying F * by S 

—x —X 

is tantamount to rescaling the rows of ( D T ) by the matrix D and then 

mm 2x 2x 

post-multiplying F by ( D T ) . The matrix (D T) is simply the matrix 

of direction numbers defining the termini of the long reference vectors 
with respect to the initial fixed axes . 

The variance covariance matrix of the second iterative stage long 
reference vectors is computed directly from Y through pre- and post- 
multiplication of Yj by S . (See Table 30 for Y*; for continuity Y ^ is 
reported in Table 31.) 
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Table 30 

Variance Covariance Matrix of the Second Iterative 
Stage Long Reference Vectors - Matrix y* 



A 2 

B 2 

°2 



A ' 

2 



B 2 



C 2 



0.522 

0.066 

0.097 



0.066 

0.511 

■0.079 



0.097 

■0.079 

0.541 



Table 31 

Intercorrelations of the Unit Length Reference Vectors 
As Determined By The Second Iterative Stage - Matrix Y t 



r* 

2 



~4x r 




Equation 41 may be further simplified algebraically as: 



*2 “ ^ 



[41] 



[42] 



Following logically from Equation 42 is the computational equation for the 
variance covariance matrix of the second iterative stage long primary 



vectors. (See Table 32.) 

z 2 = s - 4 



[43] 



From the discussion of the first two iterative stages of the 
simplified obliquimax it may be inferred that through an application of 
the laws of exponents any iterative stage can be developed directly from 
the first iterative stage without computing intermediate iterative stages. 
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Table 32 



Variance Covariance Matrix of the Second Iterative 
Stage Long Primary Vectors - Matrix Z* 

2 





A primary 

-j 


B ’^-primary 

C* 


C primary 

u 


A ’ -primary 

u 


2.031 


-0.327 


-0.411 


B’ -primary 

u 


- 0.327 


2 057 


0.359 


C ’^-primary 


-0.411 


• 0.359 


1.975 



Assume that one wishes to determine the long reference structure matrix 

associated with the U-th iterative stage. 

F* r FT(T'D~ X T) „ (T’D~ X T) (: T’D~ X T ) . . . (T'D~ X T)„ 
u j 2 3 u 

F* = FT(T'D~ X T) U = FTS U 

In? 

F* = FT(T , D~ UX T) = FTS U 

In * 

„-ux m __ „u 
p = FD T - FTS 

u 



[44] 



Because the matrix of direction numbers, S f is symmetric and because 

it is a constant at all iterative stages the matrix (FT) need only be 

u 

post-multiplied by (S ) to determine the w-th iterative stage long 



reference vectors. For the m-th iterative stage the matrix F 



(m-1) 



is 



simply post-multiplied by S to determine j*. 

The variance covariance matrix of the w-th iterative stage long 

reference vectors may be determined in a similar fashion. Assuming that 

the variance covariance matrix is computed by a pre- and post-multiplication 

of the previous iterative stage variance covariance matrix by S : 

Y* - (T’D~ X T ) n .. . (T'D~ X T) (T’D~ X T) 0 T'D~ X T) (T'D~ X T ) . . . (T'D~ X T) ; 
u u 2 112 u 



* - X u , -X . u 

Y - (T'D T) (! VD T) ; 

Y* u = Cz u d~ x t) 2u = s 2u ; 

Y* = (T'D~ 2uX T) = S 2U . 
u 



[45] 
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Thus, the variance covariance matrix of the u ~ th iterative stage long 
reference vectors is just the matrix of direction numbers raised to the 



2u - th exponential power. Following from Equations 42, 43 and 45 the 
variance covariance matrix of the u - th iterative stage long primary 
vectors is: 

z * « (t'd 2 ux t) = s~ 2u . 

u 



[46] 



2 

Using Equations 24 and 27 the diagonal rescaling matrices D ' ^ and 
2 ^ ^ 

D may be determined from I and Z respectively. The columns of matrix 
u2 u u 

* „ 2 

F as defined by Equation 44 may be rescaled by D n to form the reference 
u ul 

structure matrix, F^, or they may be rescaled by D % to form the primary 
pattern matrix, W ’ 



u 



Table 33 



Diagonal Rescaling Matrix Necessary to 
Rescale F * to F - Matrix DZ~ 

a u . al 



A f 

A 2 


B 2 


C 2 


1.383 


.000 • 


.000 


.000 


1.399 


.000 


.000 


.000 


1.360 



Table 34 



Diagonal Rescaling Matrix .Necessary 
To Rescale F to - Matrix F 



A ^-primary 


B ^-primary 


CL-primary 


1.43 


0.00 


0.00 


0.00 


1.44 


0.00 


0.00 


0.00 


1.41 
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-2 -2 

The matrices Z? and Z? are reported in Tables 33 and 34 

21 22 

respectively. The two column rescalings of F to form F and W are 

u 6 u 

reported in Tables 29 and 35 respectively. The matrices Y and Z have 

Li Li 

-1 - 1 

been rescaled by D ^ and D ^ respectively to form Y and Z ^ which are 
reported in Tables 31 and 36 respectively. 

Table 35 

Primary Pattern Matrix As Determined By 
Second Iterative Stage - Matrix W 

Li 





A -primary 

Ci 


B -primary 

u 


C^-primari 


1 


-0.027 


-1.001 


-0.028 


2 


0.044 


-0.055 


-0.976 


3 


0.999 


0.027 


-0.007 


4 


0.019 


-0.514 


-0.767 


5 


0.869 


-0.369 


0.000 


6 


0.842 


0.004 


-0.393 


7 


-0.021 


-0.756 


-0.550 


8 


0.379 


-0.864 


-0.031 


9 


0.493 


0.003 


-0.781 


10 


0.001 


-0.655 


-0.662 


11 


0.648 


-0.667 


-0.013 


12 


0.679 


0.010 


-0.622 


13 


-0.036 • 


-0.988 


-0.059 


14 


0.104 


0.001 


-0.955 


15 


0.972 


0.006 


0.037 


16 


0.741 


-0.276 


-0.369 


17 


0.287 


-0.670 


-0.489 


18 


-0.016 


-0.976 


-0.010 


19 


-0.006 


-0.105 


-0.943 


20 


0.987 


0.029 


-0.017 



Table 36 

Intercorrelations of Unit Length Primary Vectors 
As Determined By The Second Iterative Stage - Matrix Z 





A -primary 

u 


B -primary 

Li 


C2~primary 


A ^-primary 


1.000 


-0.159 


- 0.206 


B ^-primary 


-0.159 


1.000 


0.177 


C -primary 

u 


-0.206 


0.177 


1.000 
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Comparing Tables 29 and 31 with Tables 13 and 14 of Thurstone’s 
solution it is evident that the second iterative stage solution as 
determined by the simplified obliquimax closely corresponds to the final 
oblique solution as determined by Thurstone. For this illustrative 
example a reasonable estimate of p was ( ux - .40). For any oblique 
solution the final estimate of p as determined by the u~ th iterative 
stage will be ux. 

In the simplified obliquimax the symmetric matrix of direction 
numbers, S , is used in an algebraically similar manner as Thurstone’s 
matrix H (Equations 6-8). Thurstone post-multiplied the m- th reference 

fftvt 

structure matrix- by to obtain the w-th reference structure matrix. 

In the simplified obliquimax the (u-1)- th long reference structure matrix, 
(u - 1 * m) , was post-multiplied by the symmetric matrix S to obtain the 
u- th long reference structure matrix. At each iterative stage Thurstone 
had to compute H^. In the simplified obliquimax the matrix S is a 
constant matrix that is computed prior to the iterative stages. If 
Thurstone had defined S ^ as a constant symmetric matrix and used it in 
place of Hjjjj. his solution would be algebraically identical to the 
simplified obliquimax. 

Additional Geometric Asvects of the Simplified Obliquimax Matrices S and 
(D~ P T) 

Noting that the matrix S remains constant across iterative stages it 
may be correctly inferred that the i- th long reference vector of the u-th 
iterative stage has angles of inclination with the r long reference 
vectors of the (u-1)- th iterative stage that are identical to the angles 
of inclination that the i- th long reference vector of the (u-1) -th 
iterative stage has with the r long reference vectors of the (u-2)- th 

« 

iterative stage. If in Figures 10, 11 and 12 the variable points were 
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eliminated and the scalings of the axes were eliminated the remaining 
vectors would be representative of the plot of the new long reference 
vectors of the previous stage (see Figures 13, 14 and 15). It is 
prudent to keep in mind that even thuugh the reference vectors are 
oblique, the plots are on orthogonal coordinate cross-section paper. 

As in Thurstone's approach (see Subsequent Iterative Stages) it is only 
for conceptual purposes that thiL plotting approach is used. 

Figure 13 

Long Reference Vectors A* and B ^ Plotted With 
Respect to Long Reference Vectors A’ - and 

Projected Onto The ^ ane ^ 
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Figure 14 



Long Reference Vectors A* and Plotted With 

Respect to Long Reference Vectors A^_y and C^y, 

Projected Onto The A r X r - Plane 
J u-1 u-1 
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In Figures 16, 17 and 18 the termini of the variable vectors are 
plotted with respect to the planes defined by the initial fixed axes. 

The unit reference vectors as determined through five successive iterative 
stages (x - .20) have also been plotted and labeled in these Figures. From 
these figures, several important geometric properties of the obliquimax may 
be inferred with respect to the matrices D , D * and T» 

The angle of inclination between the long reference vectors and the 
previous iterative stage reference vectors remains constant throughout the 
iterative stages. Although it is not clearly apparent in Figures 16, 17 
and 18 there is a very systematic divergence of the reference vectors from 
the varimax axes toward the initial fixed axes from one iterative stage to 
the next iterative stage. If Figures 16, 17 and 18 were combined and 
plotted in a 3-dimensional space it would be clearly seen that the termini 
of the reference vectors follow distinct, predictable "paths" in space 
from one iterative stage to the next. Each reference vector will appear 

to converge toward the initial fixed axis associated with the smallest 

2 

value in the matrix D . Within the framework of planar plots it may be 
observed that both reference vectors will tend to diverge from each other 
but in so doing they will converge toward that fixed axis , of the two 

defining the plane, that is associated with the smaller value in the 

2 

matrix D . 

Within the framework of the general obliquimax and, hence, the 
simplified obliquimax there are well defined paths in n-dimensional space 
that the termini of the reference vectors will follow through successive 
iterations. These paths are determined a priori by the orthonormal 
matrix T. An orthonormal matrix computed from F by the quartimax or 
equamax (assume r > 2) criteria will not define paths that are the same 
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Figure 16 

Five Successive Approximations to Long 
Reference Vectors A' and B r Projected Onto The A B 
Plane Defined By The Original Fixed Axes A Q and B Q 
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Figure 17 



Five Successive Approximations to Long 

Reference Vectors A' and C ' Projected Onto The A C 

Plane Defined By The Original Fixed Axes A„ and C ° 

O o 
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Figure 18 



Five Successive Approximations to Long 
Reference Vectors B 1 and C 9 Projected Onto The B C 
Plane Defined By The Original Fixed Axes B 0 and °C° 
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as the paths defined by the varimax criterion. Aside from the somewhat 

weak algebraic rationale for the varimax transformation matrix, two years 

of empirical investigation has firmly established the superiority of the 

varimax T in the general obliquimax transformation (The varimax T has 

never failed to give good results for a properly factored set of data.). 

Although these paths are not fully understood it is known that they meet 

at the origin. At the origin all of the reference vectors will coalesce 
« 

into a point. As the termini of the reference vectors traverse these paths 
toward the origin the reference vectors become progressively more correlated. 
These paths may be discussed superficially as a group using the simple 
analogue of (rate x time = distance). The rate at which the reference 
vector termini traverse their paths as a group may be thought of as x. A 
unit of time may be thought of as an iterative stage. The total distance 
traversed is then a function of the value of x in the simplified obliquimax 
equations and the number of iterative stages. 

It may be assumed that at some relative distance p along these paths 
there are r points defining the termini of r reference vectors such that 
the number of vanishing perpendicular projections of the variable vectors 
onto these p reference vectors is a maximum. 

Assume that (p a .40), as it does in the illustrative problem. A 
choice of ( x = .01) would denote a relatively slow rate requiring more 
time to arrive at distance (p - .40) than a rate of ( x = .20). When 

% 

(x - .01) the time required to arrive at p would be ( u = 40) the number 
of iterative stages. When ( x - .20) the time required to arrive at p 
would be (u = 2) the number of iterative stages. 

Note that in Figures 16, 17 and 18 the reference vectors passed 
through p and that the number of vanishing projections onto the reference 
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vectors determined by the fifth iterative stage is less than the number 
of vanishing projections onto the reference vectors determined by the- 
second iterative stage. The reference vectors associated with (x = .25) 
may be inferred from Figures 16, 17 and 18. For (x = .25) the reference 
vectors will always, for any particular iterative stage, diverge farther 
than those reference vectors associated with ( x - .20). The first 
iterative stage for ( x = .25) would be a better approximation to (p = .40) 
than (x = .20), but on the second iterative stage those reference vectors 
associated with ( x = .20) would be the reference vectors associated with 
(p = .40) while the second iterative stage reference vectors associated 
wifh ( x = .25) would have diverged beyond the ideal reference vectors 
associated with (p = .40). Thus, it is most prudent that the estimate of 
x be smalt . (A reasonable estimate for most sets of data would be .10.) 
The larger the estimate the more probable is the possibility of jumping 
past p in a -single iterative stage. A value of ( x = .20) was chosen for 
the illustrative example as the ideal value of p was known to be 
approximately (.40) and a value of ( x = .20) would only necessitate two 
iterative stages. 

When to Stop Iterating 

Having discussed the simplified obliquimax computationally and 
geometrically it is now necessary to briefly discuss the concept of 
vanishing projections within the obliquimax framework. Ideally one would 
like some analytic criterion to assess an oblique solution, however the 
discussion of such a criterion will not be presented in this paper. (Two 
analytic criterion have been established for use with the general 
obliquimax and an additional four are being investigated.) 

Traditionally a subjective evaluation of the unit reference structure 
matrix loadings is used to judge the suitability of an oblique solution. 
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If the small entries of the u-th iterative stage reference structure 
matrix are larger than their corresponding values in the (w-l)-th 
iterative stage reference structure matrix it may be assumed that for 
the set of possible solutions that might be defined using a specific 
value of x the (w-l)-th solution is the best in terms of the number of 
vanishing projections onto the unit length reference vectors. 

Table 37 

Long Reference Structure Matrix As Determined 
By Third Iterative Stage - Matrix ]?* 

s 





A 1 

A Z 


H 


n t 

C 3 


1 


-0.038 


-0.585 


0.006 


2 


-0.004 


-0.008 


-0.588 


3 


0.589 


0.037 


0.024 


4 


-0.022 


-0.281 


-0.451 


5 


0.506 


-0.197 


0.034 


6 


0.486 


0.030 


-0.214 


7 


-0.045 


-0.431 


-0.315 


8 


0.205 


-0.497 


0.013 


9 


0.269 


0.031 


-0.458 


10 


-0.032 


-0.367 


-0.385 


11 


0.369 


-0.376 


0.027 


12 


0.383 


0.035 


-0.357 


13 


-0.043 


-0.577 


-0.014 


14 


0.034 


0.024 


-0.574 


15 


• 0.577 


0.023 


0.050 


16 


0.421 


-0.137 


-0.196 


17 


0.142 


-0.373 


-0.271 


18 


-0.029 


-0.571 


0.029 


19 


-0.033 


-0.038 


-0.568 


20 


0.584 


0.038 


0.018 



Within the framework of the simplified obliquimax the long reference 
structure matrix may be evaluated in place of the unit length reference 
structure matrix. The number of vanishing projections onto the long 
reference vectors will also be a maximum for the iterative stage 
associated with the best estimate of p. In Tables 23, 28 and 37 the 
first three iterative stage long reference vector structure matrices are 
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reported. In comparing Tables 23 and 28 the low loadings associated with 
Table 28 are smaller in magnitude than their corresponding loadings in 
Table 23. In comparing Tables 28 and 37 it may be observed that the low 
loadings in column three of Table 37 are smaller in magnitude than their 
corresponding loadings in column three of Table 28, however the low 
loadings in column one of Table 37 are larger in magnitude than their 
corresponding loadings in column one of Table 28. This phenomena of 
reduction and inflation of loadings in Table 37 may be thought of as 
"dumping” and actually represents the first stages of the coalescing of 
factors one and three. Inasmuch as all small loadings did not decrease 
in the third iterative stage long reference structure matrix it may be 
assumed that the long reference structure matrix determined by the second 
iterative stage is the best long reference structure matrix of the set of 
possible long reference structure matrices that might be determined for 
(x ~ .20). That is to say, the third and all subsequent iterative 
stages for (x = .20) will locate the termini of the reference vectors 
beyond the ideal points defined by p. (See Figures 16, 17 and 18.) 

Summary of Section III 

In this section the simplified obliquimax was developed and explained 
through the use of an illustrative example. A symmetric matrix of 
direction numbers, S, was used in the computation of the solution matrices 
at each iterative stage. The simplified obliquimax retained the metric of 
the original factor solution throughout all iterative stages and it was 
not until the final iterative stage solution had been obtained that the 
metric was converted to that of unit length reference vectors, thereby 
providing a traditional oblique solution. 
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The essence of the simplified obliquimax is the use of the symmetric 
matrix S and the retention of the metric of the original factor solution- 
The equations presented in this section are summarized below for the 

computation of the w-th iterative stage. 

o 

1. Define D as: 

D ' 2 = diagonal [ F'F ]. [9] 

2. Determine the orthonormal T from F using Kaiser's (1958) normal 
varimax procedure. 

3. Select some small value for x , say: 
x = .10. 



4. Define S as: 



s = T’D~ X T. 



[35] 



5. The long reference structure matrix for the w-th iterative stage 
is defined as: 

F* = FT(.T'D~ 1: T) 1 (.T'D~ X T) 2 ...(T'D~ X T) u . 

Therefore: 



F* = FTS U . 

6. The variance covariance matrix for the w-th iterative stage 
reference vectors is defined as: 



[42] 






■,-x, 



-x 



-x r 



(T 'D ~T) . . (T 'D *T) (T 'D ~T) , [T 'D ~T) IT ' D~ X T ) . . . (T ; 

6 1 19 u 



- X . 



y u - 5 



2m 



[43] 



7. The variance covariance matrix for the w-th iterative stage 
primary vectors is defined as : 

Z* = (T ’D*T)- ...(T f D X T) (T ’D X T) (T ’D X T) (T • D X T ) . . . {T , D X T) ; 

M 2 1 1 2 u 

4 = * • [44] 



To convert to a traditional oblique solution it is only necessary to 
% % % 

rescale F ui X u and Z u to the metric of unit vectors. 



8. To determine the unit reference structure matrix define: 
2 • 2u 

D ui = diagonal [5 ] = diagonal . [J M *] ; 

F = F*D~ 2 
u u ul m 



o 

ERIC 



78 



I 



I 



t 



78 



9. To determine the unit primary pattern matrix define: 

' D u 2 ~ diagonal [S~ 2u ] - diagonal [3*]; 

W u = F u D ul- 

7 7 . jr 

10. Using D 2 > D u 2 » and the reference vector and primary 
vector intercorrelation matrices may be computed as: 

Y = D~ 1 S^D~ 1 — D~^Y*D~^ • 

Y U D ul S ul ul Y u D ul> 

z = v' 1 s~ 2 u d ' 1 = zrVzf;. 

u u2 u2 u2 u u2 

All solution matrices within the framework utilize the symmetric 
matrix S . The definition of S eliminates the possibility of transforming 
to singularity and the necessity of planar plots. 



Summary 

This paper was arranged into three sections; the first section being 
primarily background; the second section being primarily theoretical; the 
third section being application and theory. 

In the first section one of Thurstone f s methods of determining 
oblique transformations was discussed within the context of his classical 
box problem. This section was presented to provide a background and to 
establish the terminology and methodology used in the subsequent sections. 

In the second section of this paper certain theoretical aspects of 
the general obliquimax were discussed to provide a basis for the 
development and understanding of the simplified obliquimax. In this 
section a cursory discussion of the general Obliquimax was provided. 

Total variance was defined with respect to the perpendicular projections 
of n variable vectors onto v long reference vectors and with respect to 
the perpendicular projections of v long primary vectors onto each other. 
The equations of the general obliquimax were discussed within the 
framework of variance and variance modification. This discussion 
necessitated an algebraic comparison of the Thurstone type reference 
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structure matrix and the Holzinger type primary pattern matrix. Finally 
it was demonstrated algebraically that when the metric of the original 
factor solution is retained in place of unit vectors the Holzinger 
pattern matrix and the Thurstone structure matrix are identical. 

In the third and final section of this paper the semi-subjective 

4 

simplified obliquimax transformation was developed and discussed within 
the context of Thurstone' s box problem. A constant, symmetric matrix of 
direction numbers was discussed algebraically and geometrically. The 
general obliquimax equations were modified and re-defined within the 
metric of the original factor solution using exponential powers of the 
symmetric matrix of direction numbers and an orthogonally transformed 
version of the initial factor loading matrix. Finally the subjective 
evaluation of the "simple structure" of a solution was discussed with 
respect to the un-rescaled reference structure matrix. The discussion 
was presented within the context of the box problem and in numerous parts 
of the section comparisons and contrasts were made with the Thurstonian 
model for determining oblique transformation solutions. 

Conclusion 

The primary objectives of this paper were pedagogical. One objective 
was to provide a reliable, semi-subjective transformation procedure that 
might be used without difficulty by beginning students in factor analysis. 
A second objective was to clarify and extend the existing knowledge of 
oblique transformations in general. A third objective was to provide a 
brief but meaningful explication of the general obliquimax. Implicit in 
these first three objectives was the fourth objective which was one of 
presenting a paper that might be profitable for both the beginning student 
and the factor analytic theoretician. 

o 

ERIC 
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To these ends the simplified obliquimax was developed having as its 
basis the classic Harris and Kaiser (1964) theory of developing oblique 
transformation solutions through the use of orthogonal transformation 
matrices. Inasmuch as the Harris and Kaiser theory does encompass the 
Thurstonian approach the Thurstonian method of determining oblique 
transformations was used to provide background information and to explain 
by analogy certain aspects of the simplified- obliquimax.* 

It may be concluded that the simplified obliquimax has fulfilled the 
first three objectives of this paper. This paper has provided: 

1. a reliable, semi-subjective transformation procedure for 
beginning students in factor analysis; 

2. a clarification and extension of the existing knowledge of 
oblique transformations; 

3. a brief explication of the general obliquimax. 



*1 would like to acknowledge the assistance that X received in this 
paper from J. Whitey. His discussions with me on the geometric aspects 
of oblique solutions proved quite valuable in the overall development of 
the total paper. 
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